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Chapter 1 

Introduction 



The first sign of quantum tunneling effect, decay of metastable state, was observed by H. Becquerel in 1886 
in radioactive uranium. Explanation of this phenomenon could be based on Louis de Broglie proposition 
that matter has properties of waves which was given in 1923. In this formalism particles penetrate potential 
barriers in the same way in which light is transmitted through absorbing medium. A quantitative result could 
not be given without Schrodinger equation and Born's probabilistic interpretation of quantum mechanical 
wavefunctions introduced in 1920s. In 1927 F. Hund studied an electron in potential with two or more 
minima separated by a barrier which is classically impenetrable. This was the first theoretical approach to 
the tunneling effect. Since then, the double well potential became a standard example of a system with 
tunneling and was considered as a toy model for more complex theories. 

In a potential with multiple minima a few lowest energies are degenerate in the domain of perturbative 
calculus. Tunneling is responsible for splitting the energies. In general, this effect cannot be studied ana- 
lytically. One can give a first approximation to the splitting using semiclassical-approximation (or WKB 
approximation) which was developed by G. Wentzel, H. Kramers and L. Brillouin in 1926. In this approach, 
one finds that the ground energy is shifted by a quantity which a nonperturbative function of action of the 
classical solution in Euclidean space. This classical solution is called an instanton. 

Validity of the instanton calculus is limited to systems with widely separated minima divided by a large 
barrier of potential. Nevertheless, it has vast applications to modern field theories. In Yang-Mill theory with 
ST/ (2) symmetry, classical fields which are constant at infinity may be considered as mappings of SU (2) group 
onto itself (see [T]). These fields may have a nontrivial topological structure and thus can be divided into 
topological sectors. It turns out that all sectors can be labeled by a Pontryagin index which takes integer 
values and is a topological invariant. In each sector there can be found a topological vacuum \n). One 
vacuum cannot be continuously deformed into another due to different topological properties. Therefore, 
a perturbation analysis about one of vacua does not take into account different sectors. The Pontryagin 
index plays an analogous role to the impenetrable potential barrier in quantum mechanics. An instanton 
in Yang-Mills theory is a classical field satisfying equations of motion in Euclidean space which connects 
two such vacua. It was discovered by A. Belavin, A. Polyakov, A. S. Schwartz and Yu. S. Tyupkin in 1975 
[2] and named a BPST instanton. Presence of the BPST instanton in the theory has serious consequences 
for the structure of vacuum in QCD. Semiclassically, trajectories which satisfy equations of motion in the 
Euclidean space form instanton liquids which are ensembles of interacting instantons (see T. Schafer and E. 
Shuryak [3]). For this reason the true QCD vacuum is a superposition of topological vacua \n) multiplied by 
phases e m6 where 9 is so called vacuum angle. 

Such system can be modeled in one dimensional quantum mechanics by a periodic potential. According 
to Bloch theorem, the energy spectrum consists of continuous bands. Each energy in a band is labeled by 
an angle which appears in energy states in the phase e m0 multiplying topological vacua. Unlike in the 
quantum mechanical case, the vacuum angle enters into the Lagrangian of Yang-Mills theory and only one 
value of 8 is admissible. No energy bands are present and there is usual mass gap between vacuum and the 
first excited state. The vacuum angle in QCD is responsible for violating CP symmetry. On the other hand, 
there is no experimental evidence for CP breaking which imposes a limit on the angle \9\ < 10~ 9 . 

The anharmonic oscillator attracted a continuous attention of physicists sine 1960s (4j[5j[6l|7]. Its double 
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CHAPTER 1. INTRODUCTION 



well version with tunneling effect was extensively studied by J. Zinn-Justin, E.B. Bogoniolny and many others 
in 1970s and 1980s ED Ell- It was discovered that there are further corrections to the WKB approximation 
which can be derived from modified Bohr-Sommerfeld quantization condition |11| . There are contributions 
from n-instanton molecules (i.e. a classical path in Euclidean spacetime which is composed of instantons 
that are close to each other) which contribute to the ground energy much less then the traditional instanton. 
Secondly, each instanton molecule contribution (including single instantons) is multiplied by a series, which 
is presumably asymptotic. Moreover, as stated recently by M. Unsal [12], interactions between instantons 
can heal non Borel summability of perturbation series for potentials with degenerate global minima. 

One has to keep in mind that instanton considerations neglect perturbative contributions to energies 
which are much larger. Secondly, there are higher order corrections to instanton contributions which become 
significant at stronger coupling. It is understandable that there is a need to verify statements concerning 
instantons and see in what regime of the coupling constant the instanton picture is valid. In quantum 
mechanics there is a very efficient method, called cut Fock space method (see J. Wosiek [13]), which we 
have at hand. It originates from the variational Tamm-Dancoff method. In [H] S. Dancoff studied the 
ground energy of fields of two relativistic nucleons. He formulated an eigenvalue problem (which is the 
time independent Schrodinger equation) and proposed a trial wavefunction which would represent collapsing 
nucleons. Energy which was obtained turned out to be higher than the energy from former adiabatic 
approximation. It meant that no collapse would take place. It shows that this variational approach needs 
a good understanding of physics of the system under consideration to propose an adequate trial function. 
In the cut Fock space method one takes basis states of Fock space \n) with n smaller than a certain cut-off 
which is supposed to be large. A price for taking so many states rather than a few trial functions is that 
the calculations have to be performed numerically. On the other hand, it is very efficient at least in three 
dimensional quantum mechanics. Indeed, convergence of energies with growing cut-off has been observed 
numerically |15| . Accuracy of this method is limited only by precision of computations and the size of the 
Hamiltonian (appearing in Schrodinger equation) which is precisely equal to the cut-off. Apart from these 
limitations, this method is exact and is then a powerful tool for testing WKB approximation. 

Both, cut Fock space method and instanton calculus start from the Hamilton operator. Let us then 
shortly discuss dimensional analysis of the Hamiltonian which will simplify notation. A Hamiltonian H may 
be given in a form 

H=— P 2 +V. (1.1) 

2m 

The potential V may be given in a form 

V = mw 2 a 2 V(X/a) (1.2) 

where V(x) is a real function of a dimcnsionlcss parameter. The function V(x) can have arbitrary shape as 
long as it is bounded from below. The parameter a is scale of the potential. The semiclassical approximation 
which will be addressed in this dissertation is valid when a is large. In fact, the same limit can be obtained 
by taking h to be small which is the classical limit. Dimensions of given operators and parameters are: 

[H] = kg cm 2 s -2 , 
[P] = kg cm s , 

H = k 9, 
[a] = cm, 

The dimensionfull operators can be then substituted with dimensionless ones. 



[X] = cm, 



[til = kg cm 2 s 
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Then the dimensionless Hamiltonian is 

H = \p 2 + V = ^P 2 + a 2 V(X/a). (1.5) 

We use the symbol* only do distinguish between scaled and non-scaled potentials. 

The idea behind introducing the scale parameter a is that for large value of a minima of the potential 
are separated by a high barrier and are far away from each other. It means that they interact weakly. It 
is common for the semiclassical approximation: there is a macroscopic scale a which is big compared to 
quantum length scale \J muj/h. On the other hand, usually a small coupling constant rather than a large 
parameter is responsible for weak interactions. For this reason we introduce the coupling constant 

9=~ 2 - (1-6) 
Let us assume that the function V{x) has a minimum at x — x /a and V{xq/o) — 0. Then 

V=^p^(X-x f + O(^g) (1.7) 

and the system is a perturbed harmonic oscillator centered at xq. We will be mainly using the coupling 
constant g whereas the scale parameter a will be used for convenience in instanton calculus. 

Plan of the thesis is the following. In chapter [2] we introduce the cut Fock space method on example 
of anharmonic oscillator with X A interaction. A typical convergence of the energy levels with growing cut- 
off is presented. Then results are compared with another numerical technique - shooting method which 
employees the standard Runge-Kutta algorithm for solving differential equations. In chapter [3] we turn to 
the double well potential. The standard instanton calculus is presented and then compared with numerical 
results obtained with the cut Fock space method. As it is well known, the WKB approximation in Euclidean 
space give relevant information only on the energy splitting and not on energies themselves because it does 
not include perturbative corrections to energies. For this reason only difference of the two lowest energies is 
compared with numerical results. In chapter [4] the cosine potential in the weak coupling limit is considered. 
It is observed that in a periodic space with periodic boundary conditions "0(0) = tp(Ka) where there are 
K minima of the potential within one period of the wavefunction, there are K energies that are split by a 
nonperturbative quantity. The splittings are obtained in the WKB approximation and with the cut Fock 
space method for K = 2,3. The case of K = oo which is closest to the Yang-Mills theory is also addressed 
by both techniques. Due to discrete translation symmetry occurring in this case, the cut Fock space method 
turns out to be very efficient. In the last chapter we study the anharmonic triple well potential. It is bound 
to have different expansions about different minima. This is why there might be no tunneling between all 
three minima. In order to have it, the potential has to be fine tuned. 



CHAPTER 1. INTRODUCTION 



Chapter 2 

Anharmonic oscillator 



The aim of this chapter is to introduce the cut Fock space method and to compare it with shooting method. 
To do this we deal with the anharmonic oscillator using both approaches. For the cut Fock space method on 
needs matrix representation of the Hamiltonian. Its relatively simple structure makes this task feasible when 
one uses eigenbasis of the occupation number operator. Eigenvalues of the truncated matrix approximate 
energies of the system. The approximation gets better as the cut-off grows. Parallelly the shooting method is 
applied. Symmetries of the system allows one to reduce the Schrodinger equation Hip — Etp to a differential 
equation on the interval [0, oo) (rather than on (—00, +00)) with initial conditions at the origin. Standard 
numerical techniques can be applied to such equation. At the end of the chapter agreement of both techniques 
is checked. The two methods are placed in different representations of the Hilbert space and are thus 
completely independent. Agreement of results implies that they are correct. 
In this chapter we will be using Hamiltonian in parametrisation 

H =)-P 2 +*-X 2 + 9 -X i + c. (2.1) 
2 2 4 

For e > the system is a harmonic oscillator with frequency oj = y/e perturbed by a quartic potential |A 4 . 
For e < it is a double well potential which will be studied in detail in following chapter. The constant c is 
introduced to keep the potential (and then also energies) positive for e < 0. 



2.1 Cut Fock space method 



This method, which is performed here after [15], makes use of the Fock space in which Hamiltonian (2.1l 
has a very simple structure. 

We will now briefly present construction of the Fock space. Let a and a) be annihilation and creation 
operators respectively, satisfying canonical commutation relations 

[a,a f ] = l, [a, a] = [a f ,a f ] = 0. 

Then the vacuum state |0) is defined to be such a normalized vector, that a |0) = is satisfied. Other basis 
vectors are constructed recursively by relation \n) — -^a) \n — 1). The Hilbert space T~L spanned by all |n)'s 
is called the Fock space. From the commutation relations it follows that a \n) = y/n \n — 1). In this paper 
we constrain our calculations to a finite dimensional space Hm — lin{|n) : n < M} called cut Fock space, M 
being the cut-off. 

In order to be able to compute how the Hamiltonian acts on the basis states one has to express it in terms 
of creation and annihilation operators using X — (a^ + a), P = ^= (a^ — a) and canonical commutation 
relations. The Hamiltonian constrained to T-Lm is a sparse matrix Hm, which elements can be easily explicitly 
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calculated: 

(m\H\n) = ((n+ \)^-^ + ^,(6n 2 + 6n + 3) + c))5 mn + — ^ — ^Vn(n - l)<$ m , n _ 2 

p(n+|)-l + e , (2.2) 

H ^ V( n + 1 )( n + 2 ) ( W>+2 V ' 

+ ^ (Vn(n - l)(n - 2)(n - 3)5 m , n _ 4 + + l)(n + 2)(n + 3)(n + 4)5 ro , n+4 ) . 

For M = oo eigenvalues of the matrix i?M are energies of the system. Crucial question is the rate of 
convergence of eigenvalues to the spectrum when M grows. It has been observed [TB] that the eigenvalues 
indeed converge. An eigenvalue A(M) converges to energy E exponentially if E belongs to discrete spectrum 
and like 1/M if E is in continuous spectrum. These conclusions are based on numerical data. Behavior of 
energies against M is shown in Figure |2.1| The value of M which is needed to obtain energies with desired 
precision highly depends on the coupling constant g. 

Let us now concentrate on the structure of energies which are presented in Figs |2.1[ |2.2| The figure 
was obtained for g — 1/98, e — —1/2 and c = 49/8. For this choice of parameters, the potential V(x) = 
c + |x 2 + |x has two global minima at x = ±7 and a maximum at x — 0. One can clearly see a change of 
behavior at the value of energy E « 6. This value is height of the potential barrier, V(0) = c = 6.125. In a 
classical system, energies E > c would correspond to states which have enough energy to propagate over the 
barrier. For E < c it cannot go through the potential barrier and stays in one minimum. Then there are two 
independent states, one in the left and one in the right minimum with equal energy. This is why energies 
smaller than c are paired. Because the system is not classical but quantum, each state with energy E < c 
is a superposition of a wavefunctions localized left and right minima. The lower energy of each pair always 
corresponds to symmetric superposition and thus is a vector composed of basis vectors with even number 
of quanta \2n) exclusively. The higher one is the antisymmetric superposition which is a vector containing 
only basis states with odd number of quanta \2n + 1). These states have parity +1 and —1 respectively. In 
position representation their wavefunctions are even and odd. Since the distance between minima is big and 
the barrier is high, there is only weak tunneling between them. In the limit of infinitely separated minima 
and an infinite potential barrier, energies would be exactly degenerate. The two corresponding wavefunctions 
could be chosen is a such way that one would be localized in the left minimum and the other in the right. 
Because it is not precisely the case, the tunneling effect mixes those states causing the true eigenstates to be 
even and odd with a slight energy splitting. Representations of eigenvectors in the configuration space are 
presented in Figure |2.2| 

We will now explain the relation between energy splittings and tunneling. In the anharmonic double 
well potential, a few lowest are paired into say E% — SEk and Ek + SEk where Ek is mean energy of the 
k'th pair and the energy splitting SEk is small. These energies correspond to even and odd states 
and \tpi respectively. In configuration space, both states are very similar up to a sign in neighborhood of 
each minimum (see Fig. 2.2 1. One can than create a state \Lf.) — -^(\ip^) + If phases are chosen 

as in Fig. |2.2| then the interference is constructive around the left minimum and destructive around right 
minimum. If one evolves the state, each component acquires a different phase: 



^_ e i(E k -5E k 



V2 



(2.3) 

'(i^l+e 2 *!^)). 



At some point in time t = to, the relative phase becomes e 2lSEkt ° = — 1 and the interference is constructive in 
the right minimum and destructive in the left one. We say that the state \Lk) tunnels into the right minimum 
in time t . When the barrier grows and widens, the energy splitting becomes smaller and tunneling is slower. 
More detailed analysis of such system, with small coupling constant, will be performed in the following 
chapter. 

Let us make one more observation on convergence of energies. As one can see from Fig. |2.1| the two 
lowest energies converge at first (for M < 30) linearly rather than exponentially as it was stated before. This 
is caused by the particular choice of basis. The eigenstates of harmonic oscillator |n) have support growing 
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with n. Therefore, for small M, the states \n) with n < M do not yet explore minima of the potential V(x). 
The first state which reaches the minimum is 1 25) . A few more states are needed to form an approximation of 
a gaussian function centered at x = 7 and from M ~ 30 convergence of the energy starts being exponential. 
For higher states, e.g. the 10th, the decrease is at first exponential, which is related to forming a wave 
with energy higher than the barrier. When M is large enough and the minima are more explored, excited 
stated in both minima can be formed. Formation of those states is reflected by linear decrease of the energy. 
Than it becomes exponential again. For yet higher states, e.g. 18th, the convergence is exponential from 
the beginning and does not become linear at any M because it converges to a energy greater that c. 



E 




Figure 2.1: Energies of Hamiltonian Hm with param- 
eters e = — 1/2, g = 1/98. A few lowest energies are 
almost degenerate. 



E 




-7 7 



Figure 2.2: Eigenstates in position representation for 
parameters e = —1/2, g = 1/98. Each function is 
shifted by energy to which it corresponds. The near 
degenerate energy pairs can be observed. Note that 
the phase is chosen to be such that the low energy 
wavefunctions cover themselves on the left side. Mul- 
tiplying every second state by —1 would reverse the 
situation. 



2.2 Shooting method 



Another popular way to determine eigenvalues of H is the shooting method, 
reads: 



^f"(x) + ~x 2 f(x) + ^f(x) = Ef(x). 



The Schrodinger equation 



(2.4) 



Equation (2.4 1 has solutions for every value of E. However, only for a discrete set {E n } there exists a 



normalizable solution f(x). This observation is the essence of the shooing method. Therefore, we will fix 
the initial values at the point x = and adjust the energy parameter in order to get a solution vanishing at 
infinity. 

Since the equation (2.4 1 remains invariant under the transformation x — > —x, function f{—x) is a solution 
whenever f{x) is one. Let fo(x) satisfy (2.4 1. Then symmetric and antisymme tric parts of fo(x), namely 
f s (x) = l{fo(x) + fo(~x)) and f a (x) 
simpler initial conditions 



h{fo( x ) — fo(— x )), are also solutions of (2.4). Moreover, they obey 



fs(0) 
f s (0) 



= /o(0); 
= 0, 



/«(0) =0; 

fL(0) =/5(o). 



(2.5) 



By the virtue of Picard theorem both functions are uniquely determined by these relations. Moreover, if /q (x) 
has finite norm then / s (x) and f a (x) also have. This shows that it is possible to consider only symmetric or 
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antisymmetric functions at a time, 
following initial conditions 



We will constrain search for energies to functions satisfying one of the 



(2.6) 



/(0) =1; f /(0) =0; 

/'(0) =0, \ /'(0) = 1. 

which leaves us with E as the only free parameter. 

The question that now arises is which values of E shall be considered to be eigenvalues of the Hamiltonian 



when equation (2.4 1 is solved numerically. For a given E one may write the differential equation in the form 

9 



f"(x) = (ex' + ?-x*-2E)f(x), 



(2.7) 



It follows, that for x such that ex 2 + |x 4 — 2E < signs of f"(x) and f(x) are different. Thus, the function / 
is "accelerated" towards the zero value and the equation is stable. However, beyond this limit the equation 
begins to be unstable, so one cannot expect any good behavior of a solution computed numerically for x 
greater than some critical value. As an example, a plot of a solution of 

V(*); 

(2.8) 



/(0) 
/'(0) 



2- 

i; 
o 




obtained numerically is presented in Figure |2.3| 
The solution is supposed to be wavefunction of the 
ground state of harmonic oscillator. The equation 



(2.8 1 loses stability at x = 1 while the solution starts 
exploding towards +oo about x value of 6.5. Similarly, 



a numerical solution of (2.4 1 taken with exact value of 



energy should converge to zero on a wide interval in 
the unstable region before it diverges. Clearly, width 
of the interval depends on precision of computations. 
We may then aim for functions which get closest to 
zero together with derivative. It can be done by find- 
ing minima of the function 



Figure 2.3: Numerically computed harmonic oscilla- 
tor ground state wavefunction. Loss of stability can 
be seen near x = 6.5. 



m(E) 



- min {|/(x)| 

0<x<K 



!/'(*)! 



/ satisfies pi?) } 
(2.9) 



starting with symmetric or antisymmetric initial con- 
ditions and varying the energy. K is the upper bound- 
ary above which the function f(x) is known to be large. 
Let us first note that in the unstable region the function f(x) may tend to zero only if its first derivative has 
opposite sign than the function itself. Once the sign of the derivative agrees with sign of the function, the 
function f(x) diverges instantaneously. Thus, the simplest way to choose K is to take the first value of x in 
the unstable region for which f(x)f'(x) > 0. This is done during computations. Ex. for the ground energy 
with e = g = 1 size of step equal to 0.01 and 16 digit precision, the equation loses stability at x ~ 0.93 while 
K = 3.15. Of course, value of K grows with precision of computations. Noteworthy, decreasing the stepsize 
is not necessary do increase K, but is needed to obtain more precise values of energy. 



2.3 Comparison of the results 

Energies obtained with both methods and their relative differences are presented in Table |2.1| M was set 
to 20. As it can be seen, the relative difference for small energies is negligible and starts to grow when we 
move to higher energies. Since M is rather small, the error presumably comes from the truncation of the 
Hamiltonian matrix. Indeed, for M = 40 we get (Ep — E s )/E s < 10~ 7 up to the 7th energy. The remaining 
difference is small and can easily be explained by numerical inaccuracy and finiteness of step size used in 
Runge-Kutta algorithm used for the shooting method. 



2.4. SUMMARY 
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Ef Es (Ef — E s )/E. 



0.620927 


0.620927 


1.34 


10" 


-7 


2.02597 


2.02597 


2.67 


10" 


-7 


3.69845 


3.69845 


6.48 


10" 


-7 


5.55758 


5.5576 


3.44 


10" 


-6 


7.56842 


7.56935 


1.22 


10" 


-4 


9.70915 


9.71146 


2.38 


10" 


-4 


11.9645 


11.9697 


4.32 


10" 


-4 



Table 2.1: Energies of the anharmonic oscillator with parameters e = g = 1 obtained by the cut Fock space 
(Ep) and shooting (E s ) methods. 

2.4 Summary 

Results of this chapter confirm that approximating energy of a system by eigenvalues of a finite matrix is 
a correct approach. In fact, it is much faster than the shooting method and thus allows to reach higher 
precision. We have shown that for a double well potential energies are nearly degenerate and explained this 
fact by the tunneling effect. It is an introduction to the following chapters where we compare the numerical 
method with analytical results obtained in the semi-classical approximation. 
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Chapter 3 

Symmetric double well potential 



The symmetric double well potential, especially anharmonic oscillator with negative quadratic term, is the 
simplest and most classical example of a system with tunneling. Still, it has some nontrivial properties of 
advanced physical systems. This is not only splitting of energies due to tunneling. Perturbative series of 
ground energy is asymptotic and non Borel summable, which is a common feature in field theory. 
In this chapter we will address the anharmonic potential 

V(x) = l(x 2 -l) 2 . (3.1) 
Recall, that the Hamiltonian is then given by 

H = \p 2 + a 2 V(X/a) = l -P 2 + -^(A 2 - a 2 ) 2 . (3.2) 
2 2 oar 

Note that when one translates the coordinate system by value a, so that zero is in the left minimum, 
X — > X — a and uses the coupling constant g = a~ 2 . Then the Hamiltonian 

H = \p 2 + \x 2 fx" + ^ (3.3) 
can be viewed as a perturbed Hamiltonian of a harmonic oscillator. 



3.1 WKB approximation 

We will now sketch the method for obtaining energies of the Hamiltonian using semiclassical approximation. 
It is done after S. Coleman |17) . A detailed analysis is performed in Appendix [A] The key to this approach 
is calculating amplitudes 

(a\e- TH \±a) (3.4) 

in the large T limit. Points x = ±a are minima of the scaled potential a 2 V(x/a). First we take the minus 
sign. One may represent the unity operator in terms of bound states of the system, I = J] \E) (E\, so that 



(a\e- TH \ 



£ (a\E) (E\ e- TH \E>) (E>\ a) = £ e~ TE (a\E) (E\ - a) (3.5) 



E,E' E 



e 



TEo (a\E ) (E \ -a)+ er TE - (a\E x ) (E, \ - a) . (3.6) 



We shall now calculate the amplitude (3.4 1 using path integrals. It can be written as 

(a\e~ TH | -a)=Af [ V[x{ T )]e- s ^ x ^ (3.7) 
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Figure 3.1: A typical shape of an instanton. Position of the instanton is at r = t\. 



where x(t) satisfies the boundary conditions x(—T/2) — —a and x(T/2) — a. The Euclidean action Se[x(t)] 
is given by the formula 



S e [x(t)] 



T/2 (x 

dr -x(t) 2 + a 2 V{x( T )/a) 

-T/2 \2 



(3.8) 



We calculate the integral using saddle point approximation. Se[x(t)] is maximized by x(t) satisfying 
SS Sx(t)^ ~ ^ anc ^ boundary conditions: 



-x(t) + aV'(x(T)/a) = 0, 
x(-T/2) = -a, 
x(T/2) = a. 



(3.9) 



Expanding aV'(x/a) about the point x — —a gives an equation with the solution x(t) = —a + ce T with an 
arbitrary constant c. It means that for large negative values of r the solution remains exponentially close 
to —a, then is suddenly jumps to neighborhood of a and approaches it again exponentially fast. For this 
reason x(t) is called an instanton. We will also use term classical solution because it is a solution of classical 



equations of motion in Euclidean space. A typical shape of an instanton is presented in Fig. 3.1 For finite 
T the classical solution crosses at r = because of symmetric boundary conditions. However, for infinite 
T the problem possesses translational symmetry. We may choose any time T\ at which x(t) passes through 
zero. Ti is called position of the instanton. 

There is a very convenient way to express Euclidean action of a instanton as a function of the potential 
only. It is possible since the energy £ = \x 2 — V(x) is conserved and is zero for the instanton solution. The 
following formula is valid in the large T limit: 

S e [x(t)] 



T-s-oo 



->5 = y dx^2a 2 V(x/a) = a 2 J dx^2V(x). 



(3.10) 



Expanding Se[x(t)] around x(t) yields 
S e [x(t)} ps + 



Sn 



dr'dT 11 



S 2 S E [x(r)} 
5x(t')5x(t"] 



Sx{t')Sx(t") 



Sx(t)=0 



dTdx(r) [ - j-j + V"(x(r)/a) ) Sx(t), 



(3.11) 



where 5x(t) — x(t) — x(t). Let us denote by {A„} eigenvalues of the operator in (3.111 
d 2 



dr 



2 + V (x(r)/a) ) x n (r) = X n x n (r), 



(3.12) 



x n (±T/2) = 0. 

Eigenfunctions x n form a complete set of functions. Expanding 5x{r) in this basis 5x(r) — ^ c n x n (r) gives 

^ oo 

S e [x(t)] ~S + - A " c «- ( 3 - 13 ) 



n,m— 
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r s 









Figure 3.2: A 3-instanton. 



One may change variables from integration over all paths x(t) to integration over coefficients {c n }. Additional 
constant coming from change of variables is absorbed by J\f. 



a e 



-TH I 



oo ., 



dCr, 



/2tt 



:e -*A»«S 



(3.14) 



n=0 v 

e _5o A/'det"5 



rfr 2 



V"(x(r)/a) 



Determinant is by definition a product of all eigenvalues of relevant operator. For large T the action does 
not change significantly if one changes position of the instanton t±. It means that the system has a zero 
mode. It is reflected by the fact that the lowest eigenvalue Ao — > when T — > oo. We have then to integrate 
out the zero mode explicitly in order to avoid a divergence in the determinant. Proper normalization of the 
zero mode yields cq = T\y/So- Then. 



da 



dc 



1 = 




dri = 



/2n J V2tt J -tii V 2?r " V 2,t 

Let det' denote the determinant with the lowest eigenvalue omitted. Then, 

d 2 



-TH i 



e -s J^TM[dr\' 
V 2n 



dr 



2 +V"{x{r)/a) 



-1/2 



(3.15) 



(3.16) 



There are other approximate classical solutions which are called multi-instantons. If individual compo- 
nents of such solutions are widely separated then they can be constructed from one-instanton solution by 
gluing instantons and antiintstantons (instanton reversed in time) in a sequence. Instantons and antiinstan- 
tons must be glued anternatively so that they form a continuous function. A plot of 3-instanton is shown in 



Fig. 3.2 Let us denote an n-instanton solution by x n (r). The Euclidean action for x n (r) is simply n times 
larger than for a single instanton: S'b[x„(t)] ss uSq. Because of the boundary conditions (3.9 1 the number 



of instantons n has to be odd. In the other case, (3.4 1 taken with the plus sign, n would have to be even. 



An n-instanton has n zero modes due to translational symmetry of each instanton separately. We shall then 
integrate them out as previously. Note that we cannot change order of instantons. 



t/2 r^r- ,T/2 

\ — dr\ l \l ~ d-T2 

-T/2 V 27T J Tl V 2lT 




i/2 



7zT 



(3.17) 



Calculation of the determinant is technical and is presented in the Appendix [A] The final formula is 

-1/2 



AM det' 



dr 2 



+ V"(x n (T)/a) 



1 



e 2 



aA 



(3.18) 
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The constant A is defined as follows. Let x(t) be the one-instanton solution: 

i£(t) = aV'(x(t)/a), 
x(—oo) = —a, 
2(00) = a. 



Then asymptotic behavior of x{t) is 

x{t) fa aA ± e^ T] 



±00, 



(3.19) 



(3.20) 



and the constant A is defined to be A — ' A+A-. Contributions from all n-instantons have now to be 
summed. One has then an arbitrary many instantons which arc separated from each other. It is known as 
the dilute instanton gas approximation. It assumes that there can be arbitrary many instantons as long as 
they are separated by a time interval much larger than the size of a single instanton. This condition is not 
taken into account in the integral (3.17). However, corrections are of order T~ x . Every classical solution has 
to begin at —a and end at a so only functions with odd number of instantons contribute: 

1 00 1 / A \ 2n+l 

^-Y L_ 

JF^L 2n+l) 



a e 



-TH\ 



e 2 



n=0 

1 • V. 
= e 2 —= sinn 

\/7T 



V 



-S ^ T 



(3.21) 



The other amplitude, (a\e TH \a) is calculated in the same way. The only difference is that now only even 
numbers of instantons contribute: 



(a\e- TH \a) 



00 

— V — 

./¥ ^ (Or, 



e 2 



cosh ( e" So -^T 



2n+l 



(3.22) 



All previous considerations were valid for any potential V(x) with minima at x — ±1. For the potential 



(3.1 ) we have 



V{x) = -{x 2 -l)\ 

S = a 2 [ dx^2V(x) = a 2 \ [ dx(l - x 2 ) = ^a 2 , 
J -l 2 J_ 1 3 



x(t) — atanh(r/2). 
It is simple to observe that A = 2. Finally, 



(3.23) 

(3.24) 
(3.25) 



(-a\e- TH \-a) 



. cosh 



2„2 2a „ 

\/7T 



a e 



-TH\ 



: sinh e a 



2 n 



(3.26) 



Using (3.5 1 we read off the two lowest energies and amplitudes of corresponding eigenstates at minima. 
Expressed with the parameter g — a~ 2 they take the form 



Eq — - — 
E x = 



,-2/3 a 



2 

— — ( 



-2/3 9 



(a\E ) = (-a\Ea) = (4tt)- 1/4 , 
{^E,) = -{-a\E l ) = {^)- x/i 



(3.27) 



3.2. COMPARISON WITH THE TAMM-DANCOFF METHOD 
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Phases were chosen such that (a\Ei) > 0. 

One has to remember that the WKB approximation does not include perturbative corrections to energies, 



which are much greater than the nonperturbative terms in (3.27). These perturbative contributions are 
identical for E and E\. It follows that the instanton calculus provides us relevant information about the 
energy difference 

AE WKB = -^e- 2 ^ (3.28) 

rather than about energies themselves. 



3.2 Comparison with the Tamm-Dancoff method 



It is instructive to compare the semiclassical approximation with a numerical method which we know that 
is essentially exact. This allows us to estimate accuracy of the WKB approximation and indicate in which 



region of g it is valid. To this end we write the Hamiltonian (3.2 1 in the form 



H = \p 2 - -X 2 



| X 4 



1 



(3.29) 



Energies of the system can be computed in the way that was introduced in the preceding chapter, by 



constructing a finite matrix that will approximate the Hamiltonian. The formula (2.2 1 for matrix elements 
of Hamiltonian may be used with substitutions 



9 



9 
2' 



(3.30) 



The constant which appears at the end of the formula (|3.29| and was not present earlier has to be added 
to diagonal elements. Then the matrix elements arc 



(m\H\n) = 



(n + |)s 



(n- \)g 



'^j y/n(n- l)5 m , n - 



v / (n+l)(ra + 2)<5 m , n 



+ 2 



(3.31) 



+ |> (y/n(n - l)(n - 2)(n - 3)<J m ,„_ 4 + y/{n + l){n + 2){n + 3)(n + 4)<J m , n+4 ) . 

In the first step it is checked for which g the 
WKB approximation result AEwkb agrees with 
numerical outcome AE num . It is shown in Figure 
3. 3 1 that AE num and AEwkb coincide for g s» 0.08. 
More precisely, for g = 0.04 the relative difference 
S = {AE WK B—AE num )/AE WK B is around 7% and 
decreases when g becomes smaller. The numerical 
procedure allows to go to much smaller values of g. 

Since the energy splitting AE num is small, high 
precision computations have to be performed. We 
succeeded to reach g — 8 • 10~ 5 which requires 3600 
digits to be taken into account in all computations. 
As we expect AE num to be of the same order as 
AEwkb, one shall consider at least a few more dig- 
its than log 10 AEwkb- This precision can be easily 
obtained eg. in Mathematica by setting precision of input to desired number of digits. To keep cut-off effects 
insignificant one has to take M large enough so that AE num does not depend significantly on it. Needed 
cut-off appears to be approximately M — 1.6/g, i.e. for g = 8 ■ 10~ 5 we took M = 20000. 




Figure 3.3: Comparison of numerical and semiclassical 
energy splitting. 
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In Fig. 3.4 we present dependence of the energy splitting AE num and AEwkb (given in (3.28)) on g. 
It is seen that both approaches are in agreement in that region. However, in order to conclude consistency 
of semi-classical approximation with cut Fock space method, one has to perform a more detailed analysis. 
More sensitive way to present our results is to plot above 5 as a function of g which is done in Fig. |3.5| From 
this it can be seen, that numerical values agree with WKB approximation as g — > 0. Still, 6 is nonzero for 
finite g. The numerical method is essentially exact and AEwkb was calculated up to a coefficient which is 
1 + 0(g). Thus our numerical data on 5(g) can be used to determine further corrections to AEwkb- The 
dashed line is a function AE num = AEwkb(1 — &9 — Pg 2 — 7<? 3 ) fitted to five points of data corresponding 
to smallest values of g. The fitted parameters are 

a = 1.47916667 ± 6.8 • 1(T 7 
P = 1.36693 ± 7.8 • 1(T 4 
7 = 4.10 ± 0.18 

The fit is consistent with the result of [8 where author proposes complete structure of the non-perturbative 
contribution. According to that paper coefficients a, /3 are obtained from calculating higher order pertur- 
bations around one instanton contribution to the energy. They are |llj 

a = — « 1.479166667 
48 



P 



6299 
4608 



1.3669705 



Coefficient 7 is not given. Coefficient a is in perfect agreement, p is away form the value given in [TT| by 
3<7. This is because of higher order terms contributions which limit the accuracy of the fit. Of course better 
approximation of a, j3, 7 and determining higher order terms would be possible if we were able to reach 
smaller values of g. 
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Figure 3.5: Relative difference between numerical and 



Figure 3.4: Numerical and theoretical energy split- theoretical energy splitting, 
ting for small parameter g. 



As we have seen, the energy splitting obtained by the WKB approximation and cut Fock space approach 
agree for g < 0.04. For smaller values of g the agreement improves. Computations were done in such high 
precision and so small values of g that further corrections were easy to determine. We have given their fittes 
values and observed that they are in perfect agreement with the results available in the literature. 



Chapter 4 

Periodic potentials 



This part of the thesis will be devoted to periodic potentials. As it will be shown in the last chapter, it is much 
easier to study tunneling when minima of potential are equivalent, i.e. Taylor expansion of the potential 
about each minimum is identical up to reflection. Periodicity of potential guarantees such equivalence. The 
most natural candidate for a periodic potential is cosine function. Such potential has infinite number of 
minima. One can takle with this problem by imposing periodic boundary conditions on the wavefunction. 
Then, the number of minima can be chosen arbitrarily. In fact, one can choose to have only one minimum and 
study contribution of instantons to the ground energy. However, no splitting of energy will take place. In the 
first section, we discuss potentials with two and three minima. In such potential multiinstanton trajectories 
can go around and get to the same minimum. We say that such trajectories have nonzero winding number. 

In second section we address the problem of a periodic potential in infinite space. Even though it is more 
complicated, the semi-classical approximation can be performed analytically. On the other hand, the cut 
Fock space method is much more challenging. 

4.1 Tunneling in periodic space 

In this section we will be interested in cosine potential in a periodic space with two and three minima. On 
one hand, there is no reason why the instanton calculus should not work in such space. Indeed, the reasoning 
is the same as in previous cases and only counting number of instantons gets more complicated. The cut 
Fock space analysis will be performed in a slightly different manner. One cannot start from eigenbasis of 
the harmonic oscillator, because there is no such system in the periodic space. On the other hand, there is 
a natural basis of plane waves which one can use to express the Hamiltonian as a matrix. 
Parametrization of the potential is as follows: 

V = g^VWgx) = (1 - cos (2n^x)) . (4.1) 

The potential has minima at x = ng~ x l 2 for tieZ. We will be interested in two cases, x G (0, 2g~ 1 / 2 ) and 
x € (0, 3(? _1 / 2 ) and will make some remarks on more general case x £ (0, Kg' 1 ! 2 ) for K minima. 

4.1.1 Instanton calculus 

The Euclidean action and one-instanton solution are independent of K. They can be calculated explicitly. 
Recall that customarily for calculations in the semiclassical approximation we use the scale of the potential 
a = g~ x l 2 rather than the coupling constant g. 

f 1 a 2 f' 1 2 

S = a 2 dx y / 2V(x) = -=- / dxy/i - cos(2ttx) = -^a 2 . (4.2) 
Jo V 2n Jo k z 
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The one-instanton solution satisfies the equation 



—x(t) + — sin (27rx(r)/a) = 0, 
2ir 

x[— oo) = 0, 
x(+oo) = a. 

It is a simplified version of the well known Sine-Gordon equation. Its solution is 

2a 



x{t) 



arctan (e r ) , 



(4.3) 



(4.4) 



up to a shift in r. The constants A±, which are defined in 3.20 

x(t) « aA±e -|r| , 

are in this case 



as 



t — > ±oo, 



A± = lim e |r| i(T)/a=- lim e |r| — 

r— f±oo 7T r— >±oo 1 -|- & 



2 

7T 



(4.5) 



(4.6) 



and A = ^JA + A_ = | . 

The number of paths of n-instanton N n depends on periodicity of the space and the two cases have to 
be addressed separately. Let us first consider the simpler case, x € (0,2a). Then there are two minima at 
zero and a. Similarly as for the double well potential, there are t wo a mplitu des to be calculated. 

Let us first calculate (a\e~ TH \0). According to formulas (3.17) and (3.181 the contribution of n-instanton 
solution integrated over zero modes is 



1 1 



(4.7) 



It does not take into account the number of topologically different paths of n-instanton solution. Each 
instanton of an n-instanton trajectory changes the minimum in which a state is localized. Therefore, only 
odd number of instantons contributes to this amplitude. Because of the periodic boundary conditions, an 
instanton starting at can go either left or right and it ends at a. There are two admissible paths for each 
instanton, so 2 2n+1 possible paths for a 2n + 1 instanton solution. 
Taking sum over n one obtains 



a e 



-TH 



10) 



e 2 



e 2 



^ oo 

.7^ 51 



22n+l 



^ (2n+l)\ 



2 n 2 2 „ 

e~^ a -^aT 



r 3/2 



2n+l 



sinh ( e ^ a 



2 4 



r 3/2 



aT 



(4. 



Conversely, only even instanton solutions contribute to (0|e TH \0). The number of 2n-instantons starting 
and ending at is given by an analogical formula, N2 n — 2 2n where n > 0. 



(0\e~ IH \0) «e" 2 



1 OO 

7?E 



i2ii 



Pk ^ (2n)l 

71=0 K ' 

: cosh ( e 



2 n 2 2 „ 
e~^ a — r-rrraT 



r 3/2 



. 2 4 



r 3/2 



aT 



(4.9) 



We now read off the energies and amplitudes. In terms of g = 1/a 2 they are 

4 



E 
E 1 



1 

2- C 
1 

2 +C 



-2/f 9. 



-2/tt 2 9 



4 



3/2 ' 



(4.10) 
(4.11) 
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(4.12) 



a 



Figure 4.1: Wavefunctions of the two lowest energy states (solid) for potential (dashed) with two minima 
for a = 10. Amplitudes in each minimum are calculated in the semiclassical approximation. Presented 
wavefunctions are gaussians in neighborhood of each minimum which is zeroth approximation in g — > 
limit. 



Let us now discuss the three minima case. Because only the topology changes and not shape of the 
potential, action of the classical trajectory So and the constant A are given by formulas (4.2) and (4.6 1. 
The remaining part is determination of numer of paths N n . Consider a triangle with vertices denoted by 
A, B, C. Each instanton moves a particle along one of the edges. Let us denote by N n (vx,V2) the number of 
n-instanton paths starting at vertex v\ and ending at vertex V2- Due to rotation and reflection symmetry, 
following equalities hold: 

<4°) = N n (A, A) = N n (B, B) = N n (C, C) 
c« = N n (A, B) = N n (B, C) = N n (C, A) 



N n (B,A) = N n (A,C) = N n {C,B) 



(4.13) 
(4.14) 



One can easily observe that = 1 and c^ 1 ' = 0. Consider now an n-instanton path which starts at A 
and ends at B. In the first step it moves either to B or C and then there are n — 1 steps left. This fact 
is represented by equation N n (A,B) — N n -i(B, B) + iV„_i(C, B). Similarity, an n-instanton trajectory 
starting and ending at A moves to B or C in the first step and then it has additional n — 1 steps to return. 
It follows that N n (A, A) — N n ^i(B, A) + A r „_i(C, A). These two equations may be expressed as 



(i) 



.,(0) 



-n-H 



cl°) = 2c 



(i) 



Inserting (4.161 to (4.151 gives 



Taking c„ = r" o ne obt ains a quadratic equation for r which solutions are r — 2 and r = —1. Then, the 



cW - 2c (1) + 



(4.15) 
(4.16) 

(4.17) 



general solution of (|4.17|) is c« = a 2 



The other series yields 



Using the initial conditions, 



-2™ 
3 



(4.18) 



#> = > + i(-l) r 



(4.19) 
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Relation N n = cffi or N n — cj 1 ^ can be used to calculate amplitudes (q2 \e TH \qi) with q\ = q2 and qi 7^ q2 
respectively. From formula (4.7 1 multiplied by c„ and summed over n one gets the final result 



,-Tifi 



1 ^ 

,a\e ~ TH \a) = (-a\e- TH \ - a) = e"^ V ^ | , 

71=0 



1 / 1 



Vi V3 



exp 



T 
2 



T 



„ _ 2 2 2aT \ 2 , 



7T 

2 2aT\ 



(4.20) 



-TJJ 



|o> = 



(a\e- TH \-a) = (-ale-™ 



|0> = 



1 OO 



/7r * — ' n! 

n=0 



1 / 1 



exp 



T „ _ 2 a 2 2aT\ 1 
2 +2e =* ^)-3 eXPl ^- e 



T 
2 



1% 



3/2 



7T 



(4.21) 



From the expansion (92 |e TH |gi) = X)b e T£ (92I-®) C^kO one can easily read off the energies: 



En 



e 

2 



-2/ttV 



r 3 / 2 . 



2 TT 3 / 2 ^ 

in terms of g — 1/a 2 and solve equations for amplitudes of the ground state \Eq 



(4.22) 



(0\E ) = (~a\E ) = (a\E ) = 



1 



(4.23) 



which are determined up to a phase. The equations for amplitudes of the wavefunction corresponding to E\ 
are 



I (OIEO | 2 = I {-a\E x ) | 2 = I (a|£x) | 2 
(01^) (£x|a) = (a|£?i> (E x \ - a) = (-a|Ei) <£7i|0> 



Multiplying three expressions on the left hand side of (4.251 gives 



I (01^) H^) (-alE,) | 2 = - 



1 

27 



(4.24) 
(4.25) 

(4.26) 



which is clearly a contradiction. Indeed, we made a wrong assumption that there is only one state of energy 
Ei while there are actually two such states, i.e. the energy E\ is degenerate. Let then \E^) and \E^) be 
the two orthogonal states corresponding to energy E\. Instead of (4.241 and (4.251 we obtain 



I (0\E[ r) ) | 2 + I (0|^ 2) ) | 2 = I (a\E[V) | 2 + I (a\E{ 2) ) | 2 

= \ { -a\E^)f + \(-a\E^)f 
(0\EW) (E[%) + (0\E[V) (E[%) = (al^ 1 ') - a) + (a\E^) (E^\ - a) 

= {-a\E^){E^\0) + {-a\E^){EP\0) = 



2 

3' 



(4.27) 



One solution to these equations is 

(-a\E{ X) ) 



(-a\E\ 



(2)\ 



1 

1 

71 



(0\E[ X) ) 

(o\e{ 2) ) 







(a|^ 1} ) 
(a|^ 2) ) 



1 

1 



(4.28) 
(4.29) 



It is not unique since there is an infinite number of possibilities in which one can choose a basis of the 
two-dimensional eigenspace corresponding to the energy E\. In our choice ) is symmetric and l-Ei 2 "*) is 
antisymmetric. 



4.1. TUNNELING IN PERIODIC SPACE 



25 




-a 



Figure 4.2: Wavefunctions of the three lowest energy states (solid) for potential with three minima (dashed) 
for a = 10. Amplitudes in each minimum are calculated in the semiclassical approximation. Presented 
wavefunctions are gaussians in neighborhood of each minimum which is zeroth approximation in g — > 
limit. 



4.1.2 Numerical computations of energy levels 

In this section we will show that energies of the Hamiltonian may be computed in a space with periodic 
boundary conditions even though we may no longer use creation and annihilation operators as in the case of 
infinite space. Construction of the matrix will be presented for general number of minima K, but computa- 
tions will be done only for K = 2, 3 for which we have derived energy levels in semiclassical approximation. 
The most natural choice of basis for finite space are plane waves. We label them as \n): 



\x\n) = —j= exp 



/ 2niriy/gx\ 
l< * V K ) 



n G 



The Hamiltonian is 



H= l -P* + V= l -P 2 +g- 1 V(JgX). 



Action of momentum operator and the potential is easy to calculate in this basis: 

P 2 \n)=g 
V\n) 



( 2im 
1 



4ir 2 g 



n) 



1 



Sir^g 



+ K) 



1 

8^ 



\n-K). 



(4.30) 
(4.31) 

(4.32) 
(4.33) 



We shall note that the Hamiltonian has a translation symmetry x — > T a x = x + a = x + 1/^/g. This is 
Zk symmetry. The basis states transform under this symmetry as follows: 



'2mri\ 
exp ( ) (x\n) 



There are K different values of exp ( 2j jp L ) so the Hilbert space can be divided into K sectors 



S k = {\k + nK),n€ Z}, 



k = 0,...K-l. 



(4.34) 



(4.35) 



The Hamiltonian can be diagonalized in each sector independently. One can see that the potential part of 
the matrix of Hamiltonian does not depend on sector Sk, i.e. 



(k + mK\V\k + nK) 



1 



2 n = m 
-1 n = m ± 1 
otherwise 



(4.36) 
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On the other hand, the kinetic part does: 

1 f 1-Kn\ 2 

(k + mK\-P 2 \k + nK)=g( — \ S nm . (4.37) 

The momentum operator depends only on the value of n 2 . Therefore, matrices of Hamiltonian are identical 
in sectors Sk and S-k = Sk-u- Energies are doubly degenerate in each sector apart from So and Sk/2- 

The Hamiltonian has also parity symmetry. However, it cannot be used together with the Z^- symmetry 
since sectors Sk are not invariant under parity transformation but transform Sk — > S-k- Thus, the only 
two parity invariant sectors are Sq and Sk/2 (unless K is odd). Parity symmetry can be used to reduce 
sizes of these two sectors into S^ and S^ (k — 0, A/2) formed by cosine and sine functions respectively. It 
computations more effective in these sectors. 

For K = 2 there are two cosine and two sine sectors. We note that the two states of lowest energies lie 
in cosine sectors Sq and S^ . This observation is based on analysis of amplitudes obtained by the WKB 
approximation and on making computations for some small g. For K = 3 there are three cosine and three 
sine sectors. There are three energies of order 1/2. One of them corresponds to function in Sq and the other 
two to functions in sectors Si and 

As before, we take a finite number of basis states and construct matrices of the Hamiltonian in each 
sector separately. Precision of computations is chosen to be such that energy splitting estimated by the 
scmiclassical approximation is seen. The cut-off for states has to be chosen experimentally. 



4.1.3 Comparison of the results 

Finally, we compare results obtained with semiclassical approximation with exact, numerical data. In Figure 
|4.3| we display energy splitting for two and three minima. As one can see, energies start to agree around 




(a) Two minima (b) Three minima 



Figure 4.3: Energy splitting obtained with WKB and numerical methods for small values of parameter g. 
Agreement of analytical and numerical approaches starts at g 0.04. 

g = 0.04 similarly as for two minima. To be more precise, in both cases relative error at g = 0.04 is around 
20%. For g — 0.011 it is already only 5%. As the energy splitting converges to very fast, it is more 
instructive to look at relative difference of energy splitting obtained with both methods. In Fig. |4.4| there 
are presented plots of (AEwkb — AE num ) / AEyy k b for cases with two and three minima. The relative 
error tends to zero. It is a proof that the semiclassical approximation and Fock space method agree. For 
small g the points align on a straight line. It suggests that there are further corrections of order g to 
AEwkb- Indeed, one can fit a polynomial to the relative difference of energy splitting. We assume that 
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Figure 4.4: Plot of relative difference of energy splitting obtained with analytical and numerical method. 
Points are values obtained from numerical computations and WKB approximation. The dashed line is a fit. 
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Figure 4.5: Wavefunctions of the two lowest energy states obtained with the cut Fock space method for 
a = 10 in the case with two minima. They are very similar to wavefunctions presented in Fig. 



4.1 



AE num = AEwkb(^ — otg — (3g 2 — 7<? 3 )- Results of fitting values a, (3, 7 are presented below. 



two minima 
a = 4.31795193±9x 1CT 7 
p = 11.2246 ±3.8 x 1(T 3 
7 = 106.4 ±3.5 



three minima 
a = 4.31795193±9x 10~ 7 
P = 11.2246 ±3.8 x 10~ 3 
7 = 106.4 ±3.5 



We have thus determined further corrections to the WKB energy splitting. However, there are no results 
in the literature to which we can compare these corrections. A very interesting fact is that the corrections 
are exactly the same. This is because perturbation about an instanton does not depend on topology of the 
space. It is an indication of the fact, that the only difference between instanton calculus for two and three 
minima, even in higher orders of g, is counting of instantons, i.e. the coefficient N n . Coefficients N n are 
purely topological. 

We have shown wavefunctions of the lowest energy states in the case of two minima obtained with numer- 
ical approach for g = 0.01 in Figure |4.5| They are in agreement with the expected shapes of wavefunctions 
obtained by semiclassical approximation. 

Note that in the case of three minima amplitudes of states achieved with WKB approximation were chosen 
symmetric or antisymmetric. In the cut Fock space approach the three lowest energies lie in sectors Sq , 1S1 
and S-i. Functions from the Sq sector are symmetric, but the two sectors S%, <5>_i contain no symmetric 
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Figure 4.6: Wavefunctions of the two lowest energy states (solid line) obtained with the cut Fock space 
method for g — 0.01 in the case with three minima. They agree with the functions obtained with the 



semiclassical method presented in Fig. 4.2 The potential is drawn with dashed line 



or antisymmetric functions. In order to obtain comparable results we have to construct wavefunctions with 
positive of negative parity. Let \ipi) be the state with energy E\ in the sector Si (which is the lowest energy 
in this sector but the first excited energy in the whole Hilbert space). The parity transformation P maps 
the sector Si into <S_i. The Hamiltonian commutes with the parity operator. Therefore, the state P\ipi) 
has the same energy as \ipi) and belongs to <S_i. We now construct two new states: 



1 

71 



(IVi> + PhM). 



\E? ) ) = -f=Wi)-P\i> 1 )) 



(4.38) 



which are both eigenstates of parity with P = +1 and P = — 1 respectively and are eigenstates of the 
Hamiltonian corresponding to the same energy E x . Phases can be chosen such that the wavefunctions are 
real. The wavefunctions are plotted in Fig. |4.6| As one can see amplitudes in each minimum agree with the 
result of WKB. 

Concluding, numerical results agree with semiclassical approximation for small g. We are able to identify 
further corrections to the instanton contributions. They are identical in both cases, potential with two and 
three minima. 



4.2 Cosine potential in infinite space 

In this section we study the cosine potential in infinite space (on K without boundary conditions). One can 
view it as K — y oo limit of cosine potential in periodic space. The potential V(x) is given by 

V(x) = t^-tzCI - cos(27ra)), iei (4.39) 

( 27 

Recall that the potential term in Hamiltonian is V = g~ 1 V( y /g~X). The feature that distinguishes this system 
from the former ones is the fact that there is infinite number of states with the energy close to E = 1/2. 
These energies form a continuous band. This is a general property of systems with infinite number of minima 
as discussed by Bloch [IS]. Also it can be seen directly in WKB calculation. In numerical computations 
we take finite K and study behavior of energies with growing K. We observe that there is an interval 
(Eq — AE/2,Eq + AE/2) with K energies. In the limit K — > oo the interval is filled - it is the continuous 
energy band. The middle of the interval Eq is perturbative value of the ground energy and Eq — >• ^ as g — > 0. 
Width of the interval is a nonperturbative function of g and vanishes for g = 0. Note that the ground energy 
of the system is E — AE/2, not E . 
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4.2.1 Semiclassical approximation 

The cosine potential has an infinite number of minima. We will be interested in obtaining all possible 
amplitudes ((to + k)a\ exp(— T H)\ma) . Because of the discrete symmetry of translation by a 

((to + k)a\e~ TH \ma) = {ka\e~ TH |0) . (4.40) 

The Euclidean action and the one— instanton solution may be obtained precisely as in the previous section, 



S = a 2 J^ dx^j^(l -cos(27ra;)) = -^a 2 , (4.41) 
x(r) = — arctan(e T ). (4.42) 

7T 



Then A = — . The number of all n-instanton paths contributing to (4.40) can be estimated in the same way 
one counts probability for the simple random walk on Z. One has to make even number of steps to move by 
an even number of minima. The number of all possible paths moving 2k sites to the right in 2n or 2n + 1 
steps is given by a Newton's binomial 

***-(„*)■ <443) 

41% = 0- (4-44) 

We modify formulas for transition amplitudes from periodic case by including proper number of paths ( 4.43[ ). 
For pT0| we find 



'/r ^ (2n)! \n+ \k\J \ 7r 3 / 2 



n=\k\ 

_t 1 / _jl„ 2 4a 
h\k\ e 



(4.45) 



e 2 -^I 9lu | e ^ a __ T 



where J a (x) is the modified Bessel function of the first kind. For odd number of minima one gets 

< fc+1) =0, (4.46) 

-1 2/,- 1.1 _ ( -" + 1 

which results in 



= „ + * + (447) 



,/ , . . rpir, . t 1 1 / 2n + l \ / 2 j 2n „\ 2,1+1 

u 71 17 V?F 4^ '.,s (2n + l)! Vn + fc + V V tt 3 / 2 J 

v n=fc (fc) v ' v 7 v 7 ( 4 4 g^ 

_t 1 / _ 2 a 2 Aa \ 

= e -3__ W| ^ e ^ _ T j, 

where fco(fc) = fc for k nonnegative and ko(k) = — k — 1 for k negative. Summing up, we have 

<fto|e- M |0> = eS±I w (e-^-^T) . (4.49) 



'7T 



It is a known fact that for integer k 

I k (x) = - [ d8e xcos{e) cos(k6). (4.50) 



7T 
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After combining (4.491 and (4.501 one gets 



{At, 



„-TH 



|o> 



£ d9 cos(ke) exp (-T Q - cob^T*"" 



or in the g — a 2 variable: 



(fee 



-TH 



|o> = 



7T 



3/2 



<20cos(fc0)exp l-T[~- cos(e)e- 21 ^ 9 '-^4 — 

2 n ' y/9 



(4.51) 



(4.52) 



The energies form a continuous band, so one should use expansion with integral over energies rather than 
sum 



(ka\e- IH \0) = dE (ka\E) (E\Q) e 



-TE 



Energies are parameterized by an angle 9: 

E{6) = - ~cos{e)e- 2 t* 2 V- 



and they form a band which has width 



AE 



ft 3/2 y/9 



-2/ir 2 g_ 



9 e (0,tt) 



ft 3/2 V9' 



(4.53) 



(4.54) 



(4.55) 



Energy states are improper (not normalizable) which is a typical situation for continuous spectrum. Their 
amplitudes at each minimum can be read off from (4.521: 



(ka\E(9)} = cos(fc0). 



(4.56) 



We know that in neighborhood of each minimum the wavefunction (x\E(9)) is approximately a gaussian of 
width 1 . We can write a general expression 



(x\E(9)) = cos{x6/a)<f>(x) 
where <fi( x ) is a periodic function with period a and 

cos(x9/a)4>(x) w e~ x2/2 for |i| < a/2. 



(4.57) 



(4.58) 



I{9/na is not a rational number then one can find such integer numbers I, I' that the equation 19 /a = 21' it— tt/2 
is satisfied with arbitrarily precision. Then 



(x + la\E(9)) = cos(x9/a + l9/a)<fi(x) « sm(x9 / a)(j){x). 
Due to the discrete translational symmetry there is another eigenstate of the Hamiltonian 

(x\E{9)} 1 = sm(xe/a)4>(x). 
Because \E(9)) and \E(9)) 1 correspond to the same energy, their superposition 



(x\E(0)), 



^±ixB I a 



tlx). 



(4.59) 
(4.60) 
(4.61) 



is also an eigenstate corresponding to the same energy. In fact, all eigenstates of the Hamiltonian - also for 
rational 9/na can be given in the form (4.611. The wavefunction (4.61) is also an eigenfunction of lattice 
translation x — > x + a corresponding to the eigenvalue e ±l8 . This is in agreement with the Bloch theorem 
which states that eigenstates of the Hamiltonian with periodic potential can be always given as lattice 
translation eigenstates. They are a plane waves e tkx modulated by a periodic function <p(x). The lattice 
momentum k is inside the Brillouin zone k £ (— 7r/a, ft/a) and energy of such state is E = Eq — A cos(#). In 
our case k = ±9 /a and A = |A_E. 

Concluding, the Hamiltonian H can be diagonalized simultaneously with the translation operator T a x — ¥ 
x + a. Energy E(9) corresponds to two points in spectrum of T a , namely e and e~ t6 . 



4.2. COSINE POTENTIAL IN INFINITE SPACE 



31 




Figure 4.7: A typical wavefunction (solid) for the cosine potential (dashed). An envelope of the gaussians is 
shown (dotted). For this function we chose 9 — 1. 



4.2.2 Comparison with the Tamm-Dancoff method 

In this section we will show how to obtain energies for the cosine potential with the cut Fock space method. 
One can do it straightforwardly as it was done in Section [3.2[ i.e. write the Hamiltonian as a matrix in basis of 
the harmonic oscillator and then diagonalize it. This approach has two major disadvantages. Firstly, matrix 
of the Hamiltonian is no longer sparse. Secondly, we do not know explicit formulas for matrix elements so 
we have to expand the cosine in power series and then sum it as far as it is possible which is very time 
consuming. For this reasons we look for another approach. The crucial point is to properly exploit the 
translational symmetry x — > x + 1/^Jg. In the infinite space the translation operator T g -i/zX — x + 1/ y/g 
has a continuous spectrum {e 1 ^ : (j> G (0, 27r)}. It does not have proper eigenspaces which could be used to 
diagonalize the Hamiltonian. A remedy to this problem is to impose periodic boundary conditions. It means 
that we have to take the same Hamiltonian as in the previous section but with large K. Matrix elements of 
the Hamiltonian in sector Sk are given by formulas: 

(k + mK\P 2 \k + nK) 
(k + mK\V\k + nK) 

Independently of the number of minima, matrix of the Hamiltonian (Hk) mn is tridiagonal in each sector 
Sk and one can use a fast algorithm for sparse matrices like Ardnoldi's iteration in order to find lowest 
eigenvalues. Moreover, if one increases the number of minima K and the cutoff simultaneously, so that the 
cutoff is constant in each sector Sk, then the time of computations grows linearly and needed amount of 
computer memory is constant! It proves that it is possible to reach high number of minima with low cost 
of computational time. Additionally, we know already that energies in sector iS^ are the same as in the 
sector S-k so it suffices to find only half of the eigenvalues. One shall remember that this approach is so 
good mostly because of the very special form of the potential. For any periodic potential which is not a 
polynomial in sines and cosines, matrix of the Hamiltonian will not be sparse. Still, using this approach one 
would always benefit from the symmetry. 

We will study two aspects of the spectrum. One is width of the lowest energy band and the other is 
distribution of energies in this band. In Fig. |4.8| we plot eigenvalues of the Hamiltonian against K. One can 
see that the lowest energy E m i n is the same for each K while the highest energy E max is constant only for 
even K. For odd values of K the highest energy is lower than for even K's but converges to the value for 
even K as K — > oo. This fact can be understood on both grounds, the semiclassical approximation and cut 
Fock space approach. 



- (k + mK\P' 2 \k + nK) + (k + mK\V\k + nK) , m, n G Z 
27r(fc + nA) x 2 



8n 2 g 



K ) (4.62) 

2 n = m 

— 1 n = m ± 1 

otherwise 
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Figure 4.8: Energies of the Hamiltonian plotted against growing K. Coupling constant is g 
lowest and highest energies for each K are marked. 
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Let us consider formulas (4.55) and (4.561. We see that the lowest energy E{6) is for 9 = 0. Amplitude 
of the wavefunction is the same in each minimum: (ka\E(0)) oc cos(fc • 0) = 1. This can be realized in space 
with periodic boundary conditions with arbitrary number of minima. In contrary the higher energy if for 
6 = ir. Then, the amplitude of the wavefunction alters: (ka\E(Tr)) oc (— l) fc . It is possible in space with even, 
but not with odd number of minima. Eg. for three minima we get: 



(E(n)\0) = - (E(n)\a) = (E(n)\2a) = - (E(n)\3a) = - (E(n)\0) 



(4.63) 



This is why for odd number of minima the highest energy in the band is lower than in the case of even 
number of minima. 

In computations for the cut Fock space method we observe that the lowest energy in sectors Sk grows 



with k for k € (0, if/2) and then decreases as k goes from K/2 to K — 1. From the explicit formula (4.621 



one can see that the matrix of the potential part does not depend neither on the sector nor on number of 
minima. The part that depends on these values is the momentum part. However, in the Sq sector 



{P 2 ) mn = (mK\P 2 \nK) = g f 27 ^!) ,) = ,( 2tt/, )" <) 



and in sector S 



K/2 



(P 2 ) mn = (K/2 + mK\P 2 \K/2 + nK) = g ( 27r (^+ nK ) ) §mn = g (27r(n + l/2) f ^ 



(4.64) 



(4.65) 



so we see that they do not depend on number of minima. This is why the lowest and highest energies in the 
band are constant for even K. For odd K there is no sector Sx/ 2 and for this reason the highest energy is 
lower than for even K. 

In order to compare width of the energy band it is enough to study energies in sectors Sq and Sk/2 which 
are independent of K as long as it is even. Therefore, they are the same as in the previously considered case 
K = 2. Semiclassical prediction of the width of energy band is 



AE 



r 3/2' 



(4.66) 
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Figure 4.9: According to the WKB approximation, dependence of the energies in the lowest energy band 
on the angle 9 yields 2(E(9) — Eq)/AE = cos(0). This plot shows convergence to the cosine function for 
small g. We put E = E(n/2). Value of E is determined numerically and AE is the result of semiclassical 
approximation. 



is also exactly the same as the energy splitting in space with two minima. Then comparison of the semi- 
classical prediction with numerical computations is presented in Fig. |4.4a| which was obtained for potential 
with two minima. 

There is yet another property of the semiclassical approximation which we can compare with numerical 



results, namely energy dependence on the angle 9. Consider the formula (4.521 for k — 



(0\e- TH \0) = -^J o d9e- TE ^ (4.67) 

Recall that the energy state \E(9)} may be decomposed into eigenstates of translation operator T a cor- 
responding to eigenvalues e and e~ 10 . Basis states of a sector Sk are also eigenstates of T a corresponding 
to eigenvalue e l8k = e 27Tlk / K . For large K points e l9k cover the unit circle uniformly. In each sector Sk there 
is precisely one energy contained in the lowest energy band. Therefore, picking some large K and finding 
the lowest energy in each sector Sk enables us to find the relation between energy and angle. Remember 
that the instanton calculus does not take into account perturbative correction, so the energy density is cen- 
tered at value Eq ^= 1/2. Now we can compare energies obtained with the cut Fock space method with 
as functions of 9 with the WKB formula E{9) = Eq — cos(9)AE /2. To do this, we plotted the energy for 
g = 0.625, 0.031, 0.016, 0.08 and K = 1000 in Fig. |4.9| One can clearly see that the energies come closer 
to the WKB (or Bloch) formula as g — > 0. 

Concluding, using the semiclassical approximation we have shown that energies which are of order E = 
1/2 form a continuous band and calculated its width. We have shown that eigenstates of each energy may 
be written as eigenstates of the lattice translation operator. All above facts are in agreement with the Bloch 
theorem. Next we confirmed width of the band numerically and showed agreement in energy dependence on 
the angle 9. 



4.3 Summary 

In this chapter we were comparing results of the semiclassical approximation with the cut Fock space method 
for periodic potential V(x) = 4^2(1 — cos(27ra)) in the weak coupling limit. We studied two cases with 
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imposed periodic boundary conditions with x € (0, Kg^ 1 / 2 ) with K = 2, 3 and a separate case ifl. 

It was shown that the instanton calculus for the cosine potential is very similar to the case of double 
well. The main difference is calculating the number of topologically different n-instanton solutions with 
fixed values at infinities. The classical one-instanton solution as well as its Euclidean action were found 
analytically. In both cases, two and three minima, we found two energies in the lowest energy regime. 
However, for K = 3 one of them is degenerate which we were able to prove within the instanton calculus. 
Explicit formulas for splitting of energies were given in (4.10) and (4.221. 



On the basis of these results and observations made in Section 2.1 we state that for any finite K, 
energies which are below the potential barrier join into groups of K energies including degeneracies. Energy 
splittings in each group are nonperturbative. However, for K > 3 similar calculations are more difficult 
because recursive relations for number of n-instanton solutions become more complex. 

In the case of infinite space, there is an energy band which also has a nonperturbative width given in 
( 4.55 1 . We also remarked that the results agree qualitatively with Bloch theorem. Let ip(x) be a wavefunction 
of energy state. One may introduce an angle 8 being a relative phase of ^(x) in subsequent minima, i.e. 
ip(x + a) — e %e ij)(x). Then energy in a band corresponding to -0 is given by the formula E = E Q — cos(#). 

Thanks to the translational Zk symmetry we were able to construct a version of Tamm-Dancoff method 
which turned out to be especially efficient for large K . Using numerical analysis we were able to establish 
agreement with energy splitting for K — 2, 3 as well as width of energy band for K = oo. We also checked 
dependence of the energy on angle 8. 



Chapter 5 

Anharmonic triple well potential 



In this chapter we will analyze a triple well potential. This is the most natural generalisation of the double 
well potential. The cut Fock space method is the same as in previous chapters and the instanton calculus 
is very similar. Because of the fact that there are three minima, one expects three energies to be almost 
degenerate. Although splitting of energies is indeed small for g — > 0, we will show that it is nonperturbative 
only in very special cases. This fact indicates that one has to be careful with the WKB approximation and 
know weather it applies to certain system or not. 

A relevant potential will be now constructed. We impose conditions on the potential so it has minima at 
x = and x = ±1 and the second derivative is equal in all minima. They are 

(5.1) 



The least order polynomial satisfying these constraints is 8th order. However, then the potential is unbounded 
from below. In order to deal with this problem we impose additional requirements which are 



V(0) 


= V(±l) 


= o, 


V'(0) -- 


= V'(±l) 


-o, 


V"(0) = 


-- V"(±l) 


= 1. 



V(±l/2) = * (5.2) 
2ir z 



This condition is inspired by the function (1 — cos(27ra:)) which satisfies both, (5.1l and (5.2 1 and is 
nonnegative. The least order polynomial satisfying all conditions is 10th order. It is positive and has three 
global minima which are x — 0, ±1. 



2 V 24 27tt 2 J V 4 9tt 2 / V 8 9tt 2 J \ 6 27tt 2 
Recall that the Hamiltonian is 



H=\P 2 + -V{^gX). (5.4) 

5.1 Cut Fock space method 

In this section we perform numerical computations to obtain energies of the system. We apply the finite 
matrix method in the same way as is was done for the double well potential. Since the potential involves 
higher powers of X, the expression for matrix elements is more complex but still amenable. 



1 / 85 512 



2 V 24 277r2 / 
8 ' 9t^ J v "'~~ '"' ' V 6 27tt 



(5.5) 



r> -' r,l2 \ TO |X 8 |n)+^-^)(m|A 10 |n). (5.6) 
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If is convenient to express each power of X as in a form 

fe 

3=1 

where N — a^a is the operator of number of quanta. For fc = 4 it is 

X4 = l + 2 N+ 2 N2 + a2 {-2 +N ) + (at)2 (2 + N ) + + 4 (at)2 ' 

Then 



(5.7) 



(5.. 



\X k \n) = qkfi{n)S m ,n +^glfc,-j(")W-j 1 



n\ , . ml 

qk,j{n)O m ■n+j\ — T- 

ml V n\ 



(5.9) 



One can easily find recursive relations for functions qkj(n) and thus find explicit expressions for arbitrary 
(but finite) number of amplitudes of the form (m\X k \n). 



We construct matrix of the Hamiltonian and compute its eigenvalues 
g G (0.001,0.1) are presented in Fig 



5.1 



The three lowest energies for 
As one can see, the ground energy E differs from the two 
higher energies Ei and E 2 by a quantity of order O(g) while the difference E 2 — Ei is nonperturbative. 
Wavefunctions corresponding to the three lowest energies are shown in Fig. 5.2 . The wavefunction ipo(x) 



corresponding to the lowest energy has support in neighborhood of the central minimum only. It means that 
a state situated in the middle minimum does not tunnel to any other minimum. The other two wavefunctions 
^1,2 (x) are approximately symmetric and antisymmetric combinations of gaussians localized in left and right 
minima. A wavefunction ip{ x ) — ^{i } i{ x ) J r4 } 2{x)) is located in the right minimum. Due to energy splitting 
it tunnels to the left minimum in a finite time. We conclude that only nonperturbative energy differences 
are responsible for the tunneling effect. 
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Figure 5.1: The three lowest energies of the Hamiltonian. It can be seen that the lowest energies is lower 
than the other two by a quantity of order 0{g) 

We may ask if the energy difference E 2 — E\ is the same as in the double well potential. This question is 
reasonable since tunneling ignores the central minimum. According to calculations for double well potential 
the splitting of energies is 

E 2 -E 1 = ^-e~ So (5.10) 

where So is the Euclidean action. As it was shown in the instanton calculus for double well potential, the 
Euclidean action Sq can be expressed as a simple integral and thus can be obtained numerically: So = 
~ j\ dxy / 2V(x) 0.4036/5- The coefficient C is not possible to be determined because there is no classical 
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Figure 5.2: Wavefunctions corresponding to the three lowest energies for g = 0.025. Wavefunctions are lifted 
by values of energies they correspond to. The lowest energy correspond to wavefunction which has support 
in the central minimum only. The other two have support in the two side minima. 



trajectory like there was in the case of double well potential. A trajectory satisfying i(±oo) = ±1 and 
equations of motion in the Euclidean space needs an infinite time to get through the middle minimum. 
For this reason it is not possible to perform a calculation for energy splitting in WKB approximation. 



Nevertheless we check if the formula (5.101 holds. 
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Figure 5.3: Energy splitting between the second and third energy in triple well potential. Points represent 
numerical data and the solid line is a fit. The discrepancy for larger values of the coupling constant g is an 
effect of higher order corrections 



Numerical data of the energy difference together with a function given in (5.10 1 is presented in Fig. 5.3 



The coefficient C was fitted to several points corresponding to lowest values of coupling constant g. Points 



are very well aligned on the curve which indicates that the formula (5.101 has correct form. Little discrepancy 
visible for larger values of g is caused by higher order corrections. 

Let us now think if there is a possibility to alter the potential in such way that there is tunneling between 
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all three minima. Let us consider potential with an additional parameter 5: 

, , 1 + 5 2 / 85 512 75\ A /31 512 155\ , 
W = — ^ + (-24 + 27^-yJ^ + U-9^ + ^) S 
55 512 135 \ 8 : (13 512 ( \ ^ 10 



8 ' 9tt» 2 7 * + I T " 27^ + M 1 s 



(5.11) 



For S — it is the previous case. Moreover, it is symmetric under the parity transformation so we expect 



tunneling between right and left minimum to be preserved. Conditions (5.1) remain unchanged apart from 



the second derivative of the potential at the central minimum which is now larger: 

V"(0) = l + 5, V"(±l) = l. (5.12) 

It follows that in the zero order perturbation expansion energy in the middle minimum is |(1 + 6) and in 
side minima it is \. On the other hand, from computations done earlier we know that for S — exact energy 
in the central minimum is lower than in the side ones. It follows that if we increase 6 then the energy Eq 
grows faster than E\ and they will eventually cross. On the other hand, according to Wigner non-crossing 
theorem [19] it cannot happen because corresponding wavefunctions t\)q(x) and ipi(x) have the same parity. 

Let us consider energies Eq(S) and E*(6) as functions of S with Eq(6) being the ground energy and i?*((5) 
the first excited energy in the symmetric (parity +1) sector. Let then ipa(5,x) and ip*(5, x) be corresponding 
wavefunctions. Let f^o,*^) stand for energy eigenfunctions in nonperturbed potential, i.e. <fri(x) = ipi{8 — 
0, x). 4>o(x) is a wavefunction localized in the central minimum while the other function, <fc*(x) is a symmetric 
combination of states localized in left and right minima. 

From what was said, if the parameter S is large enough then the energy in the central minimum is larger 
than in left and right minima. Therefore, for S large enough ipo(8, x) ~ (f)*(x), i.e. the wavefunction corre- 
sponding to the lowest energy is localized in left and right minimum. The first excited energy corresponds 
to a state in the middle minimum so t/j*(5,x) ~ 4>o(x). We infer that ^q(S,x) ~ a(6)<f>o(x) + f3* (8)<f>* (x) and 
ijj^djx) rs f3(S)cf>o{x) — a(6)*4>*{x) where \a\ 2 + |/3| 2 = 1, a(0) — 1 and f3(S) = 1 for 6 large enough. Because 
the transition is smooth, there is some point S c at which |a(5 c )| = |/3(<5 C )| = -^=. It turns out that both 

coefficients have the same phase, i.e. we can choose a(S c ) — j3(S c ) = Then a wavefunction localized in 
the middle minimum, i.e. described by function (j>o(x) ^(i^o(S c ,x) +ip*(8c,x)) will eventually tunnel to 
left and right minimum. We sec that there is tunneling between all three minima. 

The question that now arises is how to find the critical value S c of the parameter 6. At S = S c both 
E and E* correspond to states which are superpositions of states in all three minima. Then the mean 
energy of a state in the middle minimum is the same as in side minima but neither of them is an eigenvalue 
of the Hamiltonian. Difference between energies Eq and E* comes from the tunneling effect only. As we 
know, the energy splitting has to be small for small g. Thus, we may obtain S c by minimizing the energy 
splitting E* — Eq. We shall remember that we neglected the lowest energy in parity —1 sector. Numerical 
computations show that it is between Eq and E*. For 8 = S c the energy E 1 * = E^ is the second excited energy 
and the first excited energy Ei is in parity —1 sector. It means that for some 8 < S c there is degeneracy of 
energy (first excited energy in even sector is equal to lowest energy in odd sector) . It is not forbidden by the 
Wigner's theorem because the energies belong to sectors with different parity. 

Minimizing the energy difference is computationally demanding and we cannot reach very small values 
of parameter g. We managed to get down to g = 0.0016. Certainly, the value of S c depends on the coupling 
constant. Plot of S c is shown in Fig. |5.4| For 6 = both energies of states in central and in side minima 
converge to g as g — > 0. It is clear that 8 C converges to 0. For smaller g one has to lift energy of the state in 
the middle minimum by smaller value in order to exceed energy of a state in left or right minimum. 

Wavefunctions are presented in Fig. |5.5| One can see that it is possible to construct a wavefunction is 
localized only in one minimum. The state \C) — ^= (|-Eo) + l-E^)) represents a state in the central minimum 

and \L) = }^(\Eq) + y/2\Ei) — I-E2)) represents a state in the left minimum. As one can see none of them is 
energy state. It follows that after certain time each of those states will tunnel to other minima. 



5.1. CUT FOCK SPACE METHOD 
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Figure 5.5: Plots of wavefunctions corresponding to the three lowest energies. The lowest wavefunctions is 
lifted by the energy it corresponds to. The other are lifted by a greater value so that they can be distinguished. 
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5.2 Instanton calculus 



The main idea of instanton calculus is very similar to cases already considered. However, there are some 
differences since V"(0) 7^ 1. We will be again interested in calculating amplitudes 



(q2\e- TH \ qi ) 



(5.13) 



with (71,(72 = 0, ±a where we put a = 1/y/g for convenience. Because of symmetry of the potential, some 
choices of q^s yield the same result. There are four different amplitudes: 



-a e 



-TH 1 



-TH\ 



-TH\ 



a) 
a) 
a) 



(0\e- T «\0) 



(5.14) 
(5.15) 
(5.16) 
(5.17) 



A single amplitude is a sum of n-instanton contributions. The number of possible topologically different 
paths of an n-instanton trajectory N n (qi 1 q2) has to be taken into account: 



(q2\e- TH \q 1 ) = J2 N n(<lu'l2){q2\e- TH \q 1 ) n 



(5.18) 



where (92 |e TH \ c ti) n is a contribution from a single n-instanton trajectory. By the same argument as for 
the double well we have 



Me' 



-TH 



.— + V"(x n (r)/a) 



(5.19) 



where x n {r) is an n-instanton solution starting at q\ and ending at q^. Calculation of the determinant is 
more onerous than for the case of double well po tential . The full reasoning is given in Appendix |A.2| Let 



us denote w| = V"(qi). The result is the formula A. 130 
fdet' -^ + V"{x n {r)/a) 



N [del ' 



d,T' 



(loiuj 2 ) 



1/4 



,1/4 



aA 



1 +LU VSq 



_ T »l+"2 

e 2 2 



(2w(l+w))™/ 4 
(2w(l+w))("- 1 )/ 4 \/r + 
(2w(l+w))(™- 1 )/ 4 % /2^ 

Note that we excluded zero modes which give an additional coefficient 



UJl 


= L02 




= = 1 




= 1, U>2 = CO 



(5.20) 



(5.21) 



Due to reflection symmetry of the potential, action of the instanton solution does not depend on starting 
and ending point of the instanton. The Euclidean classical action is 



S = a 2 / dxy/2V(x) = a 2 dx^2V{x) 



(5.22) 



This integral cannot be calculated analytically, but it is easy to do it numerically with arbitrarily high 
precision. One has to remember that the potential depends on parameter 6 and the action has to be 
computed for each value of S separately. The highest value of 5 for our range of coupling constant g is 
5 = 0.214. Computations yield 



S (S = 0)/a 2 = 0.2019, 



S (S = 0.214)/a 2 = 0.1987. 



(5.23) 



As one can see, there is only small dependence of So on 5. On the other hand, it is multiplied by a large 
number a 2 and exponentiated. For parameters in our computations including the dependence of Sq on 5 
gives a correction from 13% for the smallest coupling up to 22% for largest coupling g. 



5.2. INSTANTON CALCULUS 
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The constant A is denned as a parameter in asymptotic form of the instanton solution. Let x{t) = az{j) 
be a one-instanton trajectory going from —a to 0. Then z(t) is a solution of the equation 



z(r) = V\z(r)), 
z(—oo) = — 1, 
z(oo) = 0. 



(5.24) 



For large r there is z(r) = V'(z(t)) w w 2 z(t). From this is follows that i(r) « ^4 + e~ WT with some constant 
A + . For large negative t expansion of the potential is different and z(t) ~ A_e T . By translation in variable 
r one can make these two constants A + and A_ equal to 

A = A y(i+u) A u/(i+ u )_ (5 25) 

The constant A cannot by determined analytically and numerical estimation has to be performed. Problem 



(5.241 cannot be solved numerically because boundary conditions are set at infinity. Thus, we impose 



boundary conditions at finite time T, 



z(r) = -V'(z(t)), 
z(-T/2) = -1, 
z(+T/2) = 0, 



(5.26) 



with large T. Apart from the neighborhood of the boundaries, the solution of (|5.26| should behave in the 
same manner as solution of ( 



5.24). We expect that functions e" r i(r) and e t z(t) approach nonzero constant 
T/2 respectively. Numerical solutions of ( 5.26| confirm this supposition outside 



values for t — > T/2 and t - 

vicinity of the boundaries (see Fig. 5.6). A± are obtained by fitting constants to plateaux visible on the 



plots. The parameter A also depends on S. Our computations yield 



e~ T z(T) 
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Figure 5.6: Behavior of the one-instanton solution near boundaries. Values of A± are determined from 
height of the plateaux. 



A{5 = 0) « 0.4284, A(S = 0.214) = 0.9119. (5.27) 

One can see that the correction coming from nonzero 5 is significant. Moreover, it accumulates with correction 
of So- 

We now consider the amplitude (0\e~ TH | — a). A classical trajectory that contributes to this expression 
has to begin at —a and end at 0. It has to consist of an odd number of instantons, say 2n + 1, first of which 
goes from —a to 0. Then there follow n two-instanton pairs. Each pair starts and ends at 0. However, the 
first instanton of each pair may go through either —a or +a. This freedom of choice gives 2™ such classical 



solutions, so N2 n +i = 2™. In formula (5.201 we put u>\ = 1, UJ2 = w and change n — > 2n + 1. In the final 



formula we also include zero modes and exponent of the classical action exp(— nSo). At last, we use the g 
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variable instead of a. 

a >-2^ ze (2n+l)! I V27r i 



n=0 



1/4 / / \ 2 ' 1+1 

x uP-/\ — 4= e ^^(2 w (l+ w ))™/ 2 ^ (5-28) 



v 2tt 2 (1 + u;) 

In this amplitude there are seen two energies: 

1 + CJ 



1/4 

T 1+u 

e 2 2 



sinh ^2 3 / 4 e- 5o ^(l + uy^ATly/^j 



E Q = 
Eo — 



4 

l+UJ 



_ 2 V* e -S ^ (l+U}) -l/4 A/V g } 

+ 2 ^ e -so^ il + ul) -m A/V g, 



(5.29) 



Let us now consider (a|e _Tff | - a). Each contributing approximate classical solution consists of even 
number of instantons, say 2n + 2. The first instanton starts at —a so it has to end at while the last 
instanton ends at a and therefore it has to start from 0. In between there are n pairs of instantons which 
start and end at 0. Again, there is a constant 2" coming from the freedom of choice whether the two- 
instantons pairs go right or left, i.e. N2 n +2 = 2™. In the formula (5.201 we put ui\ = ut2 = 1 and change 
n -> 2n + 2. 



(a\e- TH \ -a) = ^= e S (cosh (V 4 e -^(1 + U y^AT/^\ lj 



(5.30) 



Now there are three energies seen: 



E = 1 --2^e- s ^(l + ^A/^, 

Ex = \, (5.31) 



E 2 = l -+2^e- s ^{l + uy^A/^g. 



For the amplitude (— a\e ™\ — a) the calculation is the same apart from the fact that there is a trivial 
(0-instanton) constant solution which gives ./V = 1. Then, 

{-a\e- TH \ -a) = J^e"* (cosh ^e^^l + ^AT/^j + lj . (5.32) 

It gives energies identical to (5.31 1. The last possibility, (0\e~ TH |0) gives again two energies, but with mean 
value tt. As we can see, various amplitudes give different values of mean of the three energies. Nevertheless, 
they all yield the same energy splitting. It shows that only the nonperturbative correction is relevant in 
semiclassical approximation in Euclidean space while the constant term is irrelevant. Still, one can read off 
amplitudes of energy states at minima (qi\Ej): 

(~a\E ) = <0|£ ) = [^-j V4 -^r-V* (a\E ) = \*-V< 

(- a \ El ) = -^tt- 1 / 4 (0|£?i) = (a\Eh) = -^f- lfi (5-33) 



(-a\E 2 ) = <0|£ 2 ) = (^) 1/4 ^^ 1/4 (a\E 2 ) = 

These amplitudes agree with plots of wavefunctions of energy states presented in Fig. |5.5| 



5.3. COMPARISON OF THE RESULTS 
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5.3 Comparison of the results 

As stated in the preceding paragraph, it is necessary to introduce a parameter 5 to rise the lowest energy 
so that the corresponding state mixes with the two higher energy states and there is tunneling between all 
three minima. Value of S was obtained numerically by finding minimum of energy splitting of the two lowest 
energies in symmetric sector. Numerical computations show that for these values of parameter S there is 
indeed tunneling between all three minima as expected. 

Basic result of the semiclassical approximation is that there is identical splitting between successive 
energies, i.e. E2 — E\ = E\ — Eq. Ratio of these two splittings is presented in Fi g. |5.7| It can be seen that 
the ratio is 1 with great precision. A plot in higher resolution is shown in Fig. |5.8| Random behavior of 
ratio of the splittings around 1 is most likely a numerical artefact. To this end we will be considering only 



1.4 
1.2 
1.0 

0.8 
0.6 



»••••• • 



0.002 



0.006 



0.01 



0.014 



Figure 5.7: Energy splitting of the second and third energy levels compared to the splitting of the two lowest 
states. As predicted by the semiclassical approximation, the splitting is the same. 
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Figure 5.8: Ratio of the two splittings in higher resolution. Deviation of (E2 — Ei)/(Ei — E$) from 1 is 
not bigger than 0.0006 and does not depend on g. Such difference can be explained by numerical errors in 
computing eigenvalues and by accuracy of determining the parameter 5. 

splitting of the two lowest energies AE = E\ — Eq. The most interesting quantity is the relative difference 
of AE obtained by both, numerical and semiclassical methods. Relevant plot is presented in Fig. |5.9| As on 
can see, the relative difference (AE num — AEwkb)/AEwkb decreases for g — > although the indication 
that if vanishes for g — is not as strong as in the cases of double well or cosine potential. This is due to 
difficulties in determining value of 5 for small couplings g. Note that unlike in previous cases this plot was 
made in linear scale because the range of the coupling constant g is limited. Definitely, the semiclassical 
approximation does not apply for g > 0.004 where higher order corrections take over. 
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Figure 5.9: Relative difference of energy splitting of the two lowest energy states. 



Summing up, we observed that for the potential with three minima, in which the minima are not equiva- 
lent, the naive supposition that tunneling takes place between all three minima fails. There is nonperturbative 
energy splitting only between the second and third energy and corresponding wavefunctions are localized in 
left and right minima. Therefore, there is tunneling only between side minima. Still, for such case WKB 
approximation can be done (which is putting 5 = in our calculations) and predicts tunneling between all 
three minima. It follows that the WKB approximation cannot distinguish whether there is a tunneling effect 
but only determine how big it is under the assumption that tunneling takes place. For nonzero parameter 5, 
when it was confirmed numerically that there is tunneling between all three minima we checked agreement of 
energy splitting obtained with both, cut Fock space and instanton methods. Agreement of both approaches 
was confirmed in the available range of data. 



Chapter 6 



Summary 



In this thesis we were dealing with one dimensional quantum mechanical systems with multiple minima. 
Such systems exhibit the well know phenomenon of tunneling. Energy states are superpositions of wavefunc- 
tions localized at each minimum. Tunneling is responsible for splitting energies corresponding to different 
superpositions. A wavefunction which is localized in one minimum is composed of several eigenfunctions 
of the Hamiltonian and has nontrivial evolution which leads to tunneling into different minima after large 
but finite time. It is well known that in the limit of weak coupling constant the WKB approximation gives 
quantitative results for energy splittings. The toy model for investigating the tunneling effect is the double 
well potential in one dimensional quantum mechanics. Despite its simplicity it cannot be solved exactly and 
only approximate solutions are available. 

In Yang-Mills part of QCD the gluon fields have minimal energy when they are pure gauge. However, 
pure gauge fields may have different topological properties and thus cannot be transformed continuously 
one into another. Therefore, each two topologically inequivalent pure gauge fields represent two different 
topological vacua. These vacua are in fact denumerable and are labeled by a Pontryagin index. This situation 
can be modeled by a periodic potential in quantum mechanics which is also considered in this thesis. 

The energy splitting in a wide variety of double well potentials was extensively studied by many authors. 
An elegant derivation for the most classical, symmetric anharmonic potential can be found in |17| . In 1980 
E. Bogomolny observed |10| that there are higher order corrections which come from instanton attractive 
interactions. After that, in 1981 J. Zinn-Justin proposed [8,9 a formula for multi-instanton contributions 
to energy. Still, it was based on observations of only few instanton interactions. With the advent of the 
increase of computer speed it became possible to verify instanton calculations numerically. Because in the 
WKB approximation it is assumed that coupling constant is small, the first question to ask is in what 
range of the coupling constant g it is still a good approximation. We have also clarified what is the rate of 
convergence of the energy splitting obtained numerically to the splitting estimated in WKB approximation. 
Another thing which was checked is how the semiclassical approximation works for other potentials with 
multiple minima. The most interesting one is a periodic potential which is the simplest possible model for 
vacuum in Yang-Mills theory. 

In the first part of the dissertation we introduced the cut Fock space method which is essentially exten- 
sion of the Tamm-Dancoff approximation. It was compared with shooting method - a classical numerical 
technique for solving eigenequations. It turned out that in the case of anharmonic potential the Fock space 
method gives satisfactory results already for small cut-offs. Moreover, it is much faster than the shooting 
method and thus is a good candidate for investigating more challenging examples. For the double well po- 
tential we presented a detailed derivation of the energy splitting in the semiclassical approximation. The 
result was compared with numerical computations for couplings spanning from 0.00008 to 1.6. To this end 
we performed computations in very high precisions ranging up to 3600 digits. For g — 0.06 the difference 
between energy splitting obtained from both methods was already smaller than 10% to become 1% for 
g = 0.006. It was shown in [5] that the instanton contribution to energies is multiplied by a perturbative 
series in g. In the WKB approximation in Euclidean space one can find only zeroth term of this series. 
Results from cut Fock space confirm that the next correction is of order O(g). Further coefficients may be 
found in calculations based on Bohr-Sommcrfeld quantization condition and several are given in |11| . Fit to 
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numerical data confirms there analytical result. 

The next step was to perform analogical analysis for periodic potential. The instanton calculus was done 
for the cosine potential on R and in periodic space where the potential had 2 or 3 minima. In the two latter 
cases there are two and three states respectively with corresponding energies split by a nonperturbative 
quantity. For the cosine potential in infinite space there is an energy band. Splitting of energies and width 
of the band obtained with Tamm-Dancoff method and in WKB approximation were compared. In all three 
cases the agreement is in accord with expectations. For g = 0.16 there is 7% deviation which decreases when 
g gets smaller. For nonzero g there are also perturbative corrections which multiply the nonperturbative 
contributions coming from instantons. We have also shown that the wavefunctions have a Bloch form. 

Finally, the anharmonic triple well potential was considered. The exact numerical result exhibited no 
tunneling from side minima to the central one, as one would expect, but only between left and right minimum. 
This is because a wavefunction located in the middle minimum has an overlap with only one energy state and 
thus its evolution is trivial. It was shown numerically that the potential may be perturbed in such way that 
there is tunneling between all three minima. However, magnitude of the perturbation has to be computed 
numerically for each g separately which is very time consuming. We performed instanton calculations for 
the perturbed potential obtaining size of energy splitting and compared with numerical data. 1% agreement 
was established for g = 0.003 and becomes better for decreasing g. 

As we have shown, the cut Fock space allowed us to investigate energies with great precisions in a wide 
range of coupling constant g. It fully confirms results of WKB approximation. The range of applicability of 
the latter was quantitatively determined. Still one has to remember that there are some nongeneric examples 
with nonequivalent minima to which the WKB approximation cannot be directly applied. 

The next step in this research would be to investigate next order instanton contributions, i.e. involving 
instanton interactions. It is known |10| that including interaction of instantons has an imaginary ambiguity 
which is of order exp(— 26*0) where Sq is action of a single instanton. On the other hand one may calculate 
perturbation expansion of the ground energy with the Rayleigh-Schrodinger theory. The series obtained 
with this method is asymptotic and Borel sum of this series again has ambiguity of order exp(— 2Sq). It 
would particularly interesting to check that these imaginary ambiguities cancel as stated by M. Unsal in [12 . 
Preliminary results are already obtained and are to be published [20 . 



Appendix A 

Instanton calculus 



A.l Double well potential 

We consider first a simple case when the potential V(x) has minima at positions x = ±1 with equal masses 
V"(±l) = 1. Recall that the Hamiltonian is 

H = X -P 2 + a 2 V{X/a). (A.l) 
Normalized eigenvectors of the Hamiltonian satisfy orthogonality and completeness relations 

(E n \E m )=5 nm , (A.2) 
Y,\E n )(E n \=I, (A.3) 

n 

where \E n ) is a bound state of the Hamiltonian H corresponding to energy E n . The potential is positive, so 
all energies are positive as well. They can be numerated in such way that Ei < E.- L+ i. The key quantity in 
the instanton calculus is the scalar product 

(a\e- TH \~a) (A.4) 
in large T limit. The operator inside the bra-ket is called an evolution operator in Euclidean time. 

(a\e-™\ -a)=J2 (a\E n ) (E n \ e - TH \E m ) (E m \ - a) = £ (a\E n ) e- TE ~ (E n \ - a) 



n,m n (A. 5) 

e~ TEo (a\Eo) (E Q \ - a) + e"™ 1 (a^) (E^ - a) 



where the last approximation is true for large T. Once we calculate (A.4 1 we simply read off energies Eq 
and Ei. 

Let e = jj^j where N is a natural number. From now on N and T will be (infinitely) large while a large 
but finite. 



(a\e- TH \-a) = (a\(e-^) N ^\-a) 

= (a\e- £H J dx N \x N ) (x N \e~ eH ...e~ eH / dxi <aci| e~ tH |- a) 

N+l 



(A.6) 



/ dxi...dx N Y[ {^i\e eH \xi_i) 
J i=i 
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where xq = —a and xn+i = a. One of amplitudes under the integral can be expressed as follows: 
(^|e- ei? |xi_i) « J dp, (xi\l - e + a 2 T/(a; 2 /a)^ (pt\xi-i) 

dpi (xi\e 



^-+a 2 V(xJa) 



\Pi) — e 



1 

2^ 



-f (Pi-i- 



1 



ipi(xi-Xi-i) 



(A.7) 



.(i(^4^i) 2 +a 2 V(^/a)) 



/2vre 



By substituting (A.7 1 to (A. 6 1 one gets 



-TH I 



N/2 



dx\ . . . dx^e 



(A.8) 



Suppose that a function ccjy(i) = a;j with £ = ^ — ~ is differentiable in the limit N — >• oo. It satisfies 
boundary conditions x(—T/2) — —a and x(T/2) = a. Then 



AT+l 

1=1 



' ^lf + a 2 V( Xi /a)) -+S e [x(t)] 



T/2 



-T/2 



dr ( ^x(t) 2 + a 2 V(x( T )/a) 



The functional Se is called Euclidean action. In the N — !• oo limit (A. 4) becomes 

(a|e- Tff |-a) = A/ 7 / D^r)^-^^ . 



(A.9) 



(A.10) 



A 7 is an ill-defined normalization constant and we will eliminate it later. Since x{r) is usually not even 
continuous, (A. 10 1 is only a formal expression. On the other hand, only those paths for which left hand side 
of (A.9 1 is small contribute to the integral and this is true when x{r) is close to a differentiable function. It 
justifies using the expression (A. 10 1. 



A. 1.1 Saddle point approximation 



We will calculate the right hand side of (A. 10 1 in the saddle point approximation. The functional derivative 
of the Euclidean action is 



(A.ll) 



The root of (A.ll) is called classical solution of the Euclidean action because it satisfies classical equations 
of motion and is denoted by x(r). To be more explicit, we write 

-x(t) +aV'(x(T)/a) = 0, 



x(-T/2) = -a, 
x(T/2) = a. 



(A.12) 



Expanding Se around x(r) gives 

Se[x(t)] = Se[x(t)} + 



dr'dT 11 



S 2 Se[x(t)} 
8x(t')5x(t"] 



Sx(t')Sx(t" 



Sx(t)=0 



where Sx(t) = x(t) — x(t). Second derivative of the action is 



S 2 S e [x(t)} 



Sx(t')Sx(t") 



= 5(t' - t" 



Sx(t)=0 



dr' 2 



+ V"(x(T')/a) 



(A. 13) 



(A. 14) 
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Note that the third derivative is 
S 3 S e [x(t)] 



8x{t')8x{t")5x{t"') 



= crHiJ - t")S(t' - T"')V (3) (x(T')/a) = 0{a- r ). 



(A. 15) 



<5x(t)=0 



Higher derivatives contain higher powers of a 1 . In the semiclassical approximation all derivatives beginning 
from the third are omitted. Then, 



(a\e- TH \-a) = M'e~ SE ^ f V[x(r)]e^ /^r)(-£+v"(s(r)/a))«»M 



The equation 



+ V"(x{r)/a) x{t) = Aar(r) 



(A.16) 



(A.17) 



together with boundary conditions x(—T/2) = x(T/2) = is a Sturm-Liouville problem and has solutions 
only for a discrete set {A n }^L of the parameter A with Ai < A^+i and limn-^ A„ = oo. Eigenfunctions 
x n {r) satisfy completeness and orthogonality relations 



T/2 



d/rXyi \t }X m (t) 8 nm: 



-T/2 



^2 X n{T)x n {T') = S(t - t'). 
n=Q 

Function Sx(t) can be represented in the basis {x n }, namely 

oo 

5x(t) = ^ c nX n (T). 



(A.18) 
(A.19) 

(A.20) 



n=0 



Finally inserting (A. 14 1 and (A.20 1 into (A. 13 1 yields 



S e [x(t)] = S e [x(t)} 

1 °° f 
+ 2 ^ c„c m / dr'dr"(5(r' - r") 

n.m— 

1 oo 

= Se[x(t)] + - Kcl + Oia- 1 ). 

n . m — 



d 2 
'd? 2 



+ V"(x(r')/a)j x n (r')x m (T") + 0(0 



(A.21) 



The expression (A. 10 1 can be now approximated by 

^-"l-^..-«.«^/ OWT , 1 .-»=.». 

OO 

= e -S E [S(r)]^"Q 



n=0 
oo 1 

-Se[S(t)]^TJ _J_ 



rfc„ e _i AnC ; 
x/2^ 



(A.22) 



n=0 



dr 2 



+ V"(x( T )/a) 



In the first step the Wiener measure V[x(t)] was replaced by a product Y[ ^f^- Additional constant coming 

from jacobian and coefficient \j\f 7 2jx implies modification of the overall normalization factor M' — > Af. The 
determinant is a notation standing for the product of all eigenvalues of an operator. 
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A. 1.2 Action of the classical solution 



We will now calculate action of the classical solution. From (A. 12 1 we have 

\^{r)) 2 = x(r)x(r) = x(r)aV'(x(r)/a) = ^-a 2 V(x(r)/a). 



Integrating both sides one obtains 



where c is a constant. Equation (A.24) transforms to 

dr = 



dx 

y/2a 2 V(x/a) + 2c 



(A.23) 



(A.24) 



(A.25) 



Boundary conditions x(±T/2) — ±a determine the constant c which turns out to be positive. Since boundary 
conditions and the potential are symmetric, there is x(0) = 0. Hence, 



T 
2 



T/2 



dx 



lo ^2a 2 V{x/a) + 2c 

The integral on the right hand side may be divided into two parts: 

r a r a r a 

I dx[2a 2 V{x/a) + 2c]~ 1/2 = dx [(x - a) 2 + 2c] ~ 1/2 + / dxr(x) 
Jo Jo Jo 



(A.26) 



(A.27) 



where r{x) is such that expressions under integrals agree. It is simple to show that r{x) < ba 1 for x <E (0, a) 
where b is some constant depending only on the shape of V(x). The first integral is calculable: 



< log 



2c 



-a + V2c + a 2 



After simple transformations one arrives to the bound 



< c < 4a 2 e 2b e- T 



for T large enough. Using (A.24 1 again, 



r-T/2 /i 

S e [x(t)] = I dr( -x(t) 2 + a 2 V(x(T)/a) 

-T/2 \ Z 
rT/2 

dr (x(t) 2 — c) 

-T/2 

dr (x(T)y/2a 2 V(x(T)/a) + 2c) -Tc 

T/2 ^ 



dx^2a 2 V(x/a) + 2c - Tc. 



(A.28) 

(A.29) 

(A.30) 
(A.31) 
(A.32) 
(A.33) 



For large T the expression Tc is negligible, so Euclidean action of the classical solution may be approximated 
as follows: 



S e [x(t)} = / dxy/2a 2 V(x/a) = a 2 / dx^2V{x) = S 



(A.34) 
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A. 1.3 Zero mode 

For large r the classical solution x(t) is close to a and for iwn the potential is approximately a 2 V(x/a) 
\{x~a) 2 . Then, 



Solution of this equation yields 



x = y / 2a 2 V(x{T)/a) + 2c « a - x(t). 



x(t) w a — Ce r . 



(A.35) 



(A.36) 



That means that for t> 1 the solution is very close to the stationary point a of the potential and contribution 
to the action is small. The same is true for r <C — 1. Thus we can say that the classical solution jumps 
from one stationary point to another in time of order 1. For this reason it is called an instanton. Major 
contribution to the action is from the region where |r| is of the order of 1 or less. It follows that a trajectory 
x Tl (t) satisfying 



x Ti (t) = x{t- ri) 



for - - |r 2 | » 1 



(A.37) 



and slightly modified near the boundaries (so that is satisfies boundary conditions) gives almost the same 
value of Euclidean action as x(t), i.e. Se[x Ti (t)] rs So- It means that the lowest eigenvalue Ao of the operator 
(A. 17 1 which corresponds to translation of the classical solution tends to zero as T — > oo. This results in a 
divergence 



dc 
\Z2n 



e 2 = oo. 



(A.38) 



For infinite T the eigenfunction corresponding to the zero mode Ao = is Xo(t) = o,4-x(t) where a is a 
normalization constant. Indeed, 



The normalization constant can be obtained from the orthogonality relation (A. 18): 

r T/2 



= 



1 = 



T/2 



cItxo(t) = a' 



T/2 



cLtx(t) = a Sq. 



(A.39) 



(A.40) 



Therefore, a — S . If x(t) = x Ti (t) is a shifted instanton then 



Sx(t) — x(t — Ti) — x(t) Ti —x(t). 

OfT 



On the other hand 8x(t) = cqXq{£). Thus, 



Co = yS Ti. 



(A.41) 



(A.42) 



Obviously t\ G (—T/2, T/2) and for finite T the integral (A.38 1 is finite since the integration limits are finite: 



dc 
V2n 



T/2 /So 



-T/2 



2tt 



-T. 



We pull the zero eigenvalue outside of the determinant defined in (A. 22 1 and arrive at 

-1/2 



(a\e ~ TH \ - a) = e- s «\ I ^TAf { det' 



~+V"(x(r)/a) 



(A.43) 



(A.44) 



The symbol det' stands for the product of all eigenvalues apart from the lowest Aq. 
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A. 1.4 Multi— instanton classical solutions 



There are other approximate classical trajectories that one has to consider. They consist of several instantons 
glued together so that the resulting function jumps from one minimum of the potential to another 2n + 1 
times. The number of glued instantons must be odd so that is starts at —a and ends at a. Let us denote 
such classical solution by X2 n +i (?). Contribution of each instanton to the action is approximately So, so 
Se[x2ti+i{t)] ~ (2n + 1)50. Let us denote by t\, t%, . . . , T2 n +i times at which x 2n+ i(r) passes zero. We call 
them positions of instantons. 

The gaussian approximation around £2,1+1 (t) yields 



— TH 1 



i) 2n+1 = e-( 2 " +1 > 5 °AAdct-5 



+ V"{x 2n+1 {T)/a) 



(A.45) 



in analogy to (A. 22 1. The operator —-§pi + V" (x2„+i (r)/a) has now 2n + 1 eigenvalues which are approxi- 
mately equal 0. They correspond to freedom of choosing positions of the instantons. We shall treat them in 
similar way as previously. They are two constraints on Tf. their order cannot be changed, i.e. Ti < Tj+i and 
they have to be separated by a distance at least 1 (which is size of a single instanton) . 



dc 
V2^ 



dc 



2n+l 



2d 



2n 



,T/2 


,T/2 


,T/2 


/ dr x 


/ dT 2 . . 


• / d,T 2n+ i 


J-T/2 


>T1 + 1 


Jr 2n + 1 



2n+l 



rp2n+ 



(2n 



(1 + 0(T^)) 



Finally, 



2n+l 



a e 



-TH I 



/2n+l 



(2n + l)\ 



-So, 



J\f det ' 



V r "(x 2n+1 (r)/a) 



-1/2 



(A.46) 
(A.47) 

(A.48) 



Here det ' stands for the product of all eigenvalues apart from the 2n + 1 lowest ones. 
A. 1.5 Determinant (part 1) 

We will now pass to calculating the determinant. To do this we recall the formula it originated from. 



Af det 



dr 2 



+ V"{x{t)/o) 



\ _1 / 2 r 

) =N' D[x(T)]e- J '^(h(Mr)) 2 + kV"(x{r)/a)x(r) 2 ) ( A49 ) 

x(t) stands for what was earlier denoted by 8x(t) and obeys boundary conditions x(±T/2) = 0. The path 
integral on the right hand side of (A. 49 1 is defined as a limit of multiple integral at time slices which can be 
transformed to scalar products of the form (xi \ . . . conversely to what was done in (A.7l. 

Af' J D[x ( T )] e - S Hl^)) 2 + lV''(*(r)/a) X (rf) = f ^ _ _ _ ^ J| ^ j - f dr( |P 2 + \V (S(r)/a)X 2 ) 



0|f(e 



-/dr(ip 2 + iy"(i(r)/a)X 2 ) 



\X = 0) 



Xi-l) 



(A.50) 



where T is time ordering operator. The time ordering appears because the last expression was defined a 
product of operators acting between time slices. The Euclidean evolution operator is defined as 



U (T2,n) = f (VA 1 dT{kPi+hV(x(r)/a)X-y 

Due to the time ordering it satisfies usual properties of an evolution operator 

U(t,t) = 1, 
U(T 2 , n ) = U(t 2 ,t 3 )U(t 3 ,ti). 



(A.51) 



(A.52) 
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As it was shown before, V"{x(t)) rs 1 unless t e (—1, 1). It indicates that U{t2,t\) may be approximated 
by the Euclidean evolution operator of the harmonic oscillator 



U (T 2l n) = t (e- A 2 M^+** a )) = f ( e -(— )(^ 2 +l^ 2 )) 
Energy states of harmonic oscillator are well known and will be used in what follows. 



(A.53) 



uiX- T 

2' 2 



£■0 



E 



e/(i,-W-i,-~ 

Z7(1,-1)C7 (-1,- 



r 



-(T/2-l)E° | R 



-(T/2-l)K,' I niO\ / tpOi tt n i\ I ciO\ /rpO\ „-(T/2- 

1)E o (E° \U (1,-1) \E°) 

e °\E )(E \U (l,-l)\E )(E \e o ______ 



E 



,-(T/2-l)E° | R 



|^}<^|c/ (i,-i)K)<k: 



I _-(T/2-l)E£, 



U 



(i? |C/o(l,-l)|^o > 

T T\ (£ °|[/(1,-1)|£ °) 



2' 2 7 <£7g| C/ (1, -1) |£7°> 
T T N 
2'~ 2 



(A.54) 



If we performed the same analysis for the classical trajectory X2n+i(i~) there would be 2n + 1 intervals of 
length 2 at which V" (x2n+i{i~)) 76 1. This results in 2n + 1 coefficients k. Therefore, 



AM det ' 



dr 2 



+ V"{x 2n +i(r)/a) 



-1/2 



: Af det 



dr 2 



-1/2 



2n+l 



(A.55) 



One may observe that Af standing in (A. 22 1 is independent of the structure of the potential or choice of the 
points a and —a standing on the left hand side of that equation. It is then justified to write 

-1/2 



Af det 



dr 2 



1 



= (x = 0\e- 1Ho \x = 0) 



(A.56) 



where Hq = \P 2 + \X 2 is the Hamiltonian of harmonic oscillator. Its ground energy and bound state are 



E 



1 
2' 



(x\E°)=n^e 



(A.57) 



Therefore, 



det 



dr 2 



-1/2 



(x = 0\e~ TH °\x = 0) = e- TE « (x = 0\E°) (E°\x = 0) + . . 

_T 1 



(A.58) 
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A. 1.6 Summing over instantons 

After summing all contributions one obtains 



(o e 



-TH\ 



^^ o (2n + l)! ^ 



V Ztt 



2n+l 



(A.59) 



There is a question of the range of summation N. Since single instantons in the trajectory X2n+i{i~) have 
to be well separated, it implies a bound for TV which is T/(2N + 1) 3> 1. On the other hand, contributions 
from multiinstantons should be included until 



(2JV+1)! I V 2tt V 



2JV+1 



< 1. 



(A.60) 



By applying the Stirling formula we find 



e ~ So \ tt-k^T -C 2N + 1. 

V Z7T 



Joining (A. 61 ) with the condition 2N the number of instantons TV has to satisfy 



'So n- 2N+1 i 
— kVAo < - < 1. 

Z7T J 



(A.61) 



(A.62) 



It is possible to choose such N since Sq oc a 2 —} oo. Terms (a|e 



— TH I 



_! with n > TV do not provide any 



information about instanton solutions. Yet they are so small that they can be added without introducing a 
significant error. Thus we set N = oo and 



{a\e~ TH \-a) = 
A. 1.7 Determinant (part 2) 





r d 2 : 




^det 




) 




~W 2 +\ 





-1/2 



sinh e 



(A.63) 



It is now sufficient to find the term kvAo- From (|A.55| it follows that 

1 



det 



Ky / \ 



dr 2 



1/2 



det'[-£^ + V"(x(r)/a)] 
Let us take a general potential W(r). Define ipx as a solution of 

LwMr) =(-J-2+ W (t)) Mr) = W\(t), M~T/2) = 0, JUa(t) 



(A.64) 



= 1. 



(A.65) 



--T/2 



For a discrete set {A„}^L of the parameter A there is ipx n (T/2) = 0. They are eigenvalues of the operator 

L\y- If W(t) — then X n = m ( r -^-) 2 ■ W(t) is a bounded function on the interval (— T/2,T/2) and can be 

regarded as a bounded operator added to the hermitian unbounded operator Thus, eigenvalues A„ of 

L w satisfy |A„— (^) 2 | < M = max(_ T / 2 ,T/2) W(t). Let us take two arbitrary potentials W^\t), W^ 2 \t) 
and define the ratio of determinants as follows: 



R(X) 



det 


d 2 

dr 2 


f Wi(t) - A 


det 


r d 2 

. dr 2 ' 


fW 2 (r)-A] 



n^ 2) -A)' 



lim 



(A.66) 
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The limit on the right hand side exists. Indeed, 



TlLoW - A) 

A) 



(i) 



JV 



nto(A 



(2) 



exp lo S 



A (1) 



A n ^ — A . 



/ N X (2) 

exp 2^ 



A (2) A 



(A.67) 



(1) 



A„ 2) | < AfW + M (2 *> so for constant A the ter m in numerator is bounded while the denominator 



behaves like n 2 . Thus, the series is convergent for N — > oo. The function i?(A) has zeros at A = A» an d 

(2) 

poles at A = A„ . Observe that for 
to 1. Therefore, 



oo each term in the product on the right hand side of (A.66) tends 
-> 1. (A.68) 



Let us define a function 



R(X) 



R'(X) 



|A|-voo 



^\T/2) 



(A.69) 



For |A| large enough the potential W(t) in (A. 65 1 is negligible compared to |A| so the function ip\(r) does 
not depend on the potential significantly. It follows that R'{X) — > 1 when |A| — > oo. The function R'(X) has 
zeros and pole in the same points as function R(X). This is merely because ip^(T/2) = <^=J> A = A« . 
Then the function 



9{x) = wm 



(A.70) 



is an entire function. It is also bounded because it converges to 1 for infinite |A|. By Liouville's theorem it 
is constant, i.e. R(X) = R'(X) for all A. 



dot 


d 2 

dr 2 




_4 1] (T/2) 


det 


r d 2 

. dr 2 " 


\-W 2 (r) 


~ 4 2) (T/2)' 



which implies 



«V Ao 



V>o(T/2) ' 



Let Vo(t) be solution of (A. 65 1 with W(r) = 1. It is easy to show that Vo(t) = sinh(r + T/2). Thus, 



The function ^)q{t) satisfies 



dr 2 



+ F"(x(r)/a)U (r) = 0, 



M-T/2) = 0, 



= 1. 



M-T/2) 

From (A. 36 1 it follows that the one-instanton solution has asymptotic behavior 
x(t) « aA±e _ ' T ' t — > ±oo. 



(A.71) 

(A.72) 

(A.73) 

(A.74) 
(A.75) 

(A.76) 



This is definition of constants A±. From (A. 39 1 it follows that the function yi{r) = -)=x{t) satisfies the 

I [ V^o 

equation (A.74 1 but not necessarily boundary conditions (A.75). Asymptotic behavior of yi(r) is 



±oo 



(A.77) 
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with B± = -r=A±. Constants B + and £?_ can be a priori different. However, we may redefine the solution 
yi (t) — > 2/i (t — t\) so that these constants change B + — » B' + = B + e Tl and £?_ — > B'_ = B_e~ Tl . We choose 
Ti in such way that B' + = B'_ = B. Wronskian W of independent solutions of a linear differential equation 
is constant. Therefore, one can normalize the other solution y 2 (r) of equation (A.74) so that is satisfies 



W = yx{r)y 2 {r) - yi(r)y 2 (r) - 2B 2 

For large |r| it reads 

2/2(r)-2/ 2 (T)=2S e H, 
Thus, the asymptotic behavior of y 2 (t) is 

y 2 (r) = ±Be |r| + C±e HT| w ±£e |T| , 



±oo. 



T — > ±00. 



Finally, ?/>o(t) satisfying both (A.74) and (A. 75) is 

implying ipo(T/2) — 1. The Green function for the differential operator £y"(2(-r)/a) is 

rfr^-f -( W )~W)2/ 2 (r) r'<r 
° IT ' TJ \ -(W)- 1 y 1 (r)y 2 (i J ) r'>r 

Therefore, the function ip\ defined by ( A.65[ ) with W(t) — V"(x(t)) fulfills the integral equation 

r T/2 



^a (t) = ^q(t) + 



T/2 



(A.78) 
(A.79) 
(A.80) 

(A.81) 
(A.82) 

(A.83) 



Function ipo(r) is added so that ijj\ satisfies initial conditions. As shown previously, Ao is close to 0. It is 
then justified to replace t/j\ with ipo on the right hand side of (A.83 1 since it introduces higher order error. 
Secondly, ip\ (T/2) = as mentioned earlier. 



« Vo(T/2) 



T/2 



dr'G(T/2,T')A Vo(r') 



-T/2 



1- Ao 



y 2 (T/2) [ T/2 , , 



4B 3 



-T/2 



dr' yi (r') (e T ^ Vl {r')+e- T / 2 y 2 (r'\ 



Ao 
4B 2 



T/2 



-T/2 



dr> (e T yi (T') 2 + yi(r')y2(r')) 



(A.84) 



1 - 



_Ao_e T 
4S 2 5 

Ap t 
4B2 6 



T/2 



T/2 



dr'x^r'f 



The term J" dT'yi{T r )y 2 (T') was omitted because it is of order 1 and is negligible compared to e T . The last 
equality is guaranteed by (A. 40 1. We obtain 



A = AB 2 e~ T . 



Finally, we get 



^X~ = V2B = J-^aA 



(A.85) 



(A.86) 



where A = ^JA+A-. 
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A. 1.8 Final result 



Inserting this result and (A. 58 1 into (A. 63 1 one obtains 

T 1 



a e 



-TH\ 



Comparison of this result with (A. 5 1 yields 



E = - - e 
2 



E 1 = ~ 2+ e 



sinh ( e- So 



So 



aA 
aA 



(A.87) 

(A.88) 
(A.89) 



A. 2 Inequivalent minima 

As demonstrated on the example of triple well anharmonic oscillator, tunneling does not take place always 
when we naively expect it. The triple well potential is 



V(x) 



1 
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..io 



(A.90) 



with the Hamiltonian H = \P 2 + a 2 V(X/a) The potential V(x) has equal derivatives V"(x) = 1 at all of 
the three minima x — 0, x — — 1, x = 1. Thus we expect that there are three energies which are close to 
E = 1/2. The next supposition is that when one locates a state in one minimum, it will tunnel to other 
minima. This statement is not true. A state with the lowest energy is localized in the middle minimum. 
Then there are two states, with even and odd parity, which are localized in left and right minima. It means 
that there is a tunneling between left and right minimum and the middle minimum is separated. Moreover, 
energy of the middle minimum is pcrturbatively different from the two higher energies, which in contrast are 
split only by a nonperturbative amount. This is clear when one notices that Taylor expansion of the potential 
V(x) about point x = is different than the expansion about x = 1. It turns out that this perturbative 
inequivalence excludes tunneling. 

As presented in chapter |5 . 1 1 when one perturbes the potential in the following way: 



V s (x) 
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(AM) 
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then the second derivative in the middle minimum is higher: V"(0) = 1 + 5 so the energy of the state localized 
in the middle minimum grows with 5. Energies of states in left and right minima grow slower because the 
parameter 5 enters expansion of V(x) about x = ±o with the term x 3 rather than with x 2 . Eventually, energy 
of the state localized in the middle minimum has to cross energies of the states in left and right minima. It 
is forbidden by the Wigner non-crossing theorem. Thus, for a particular choice of <5 the tunneling between 
middle and left and right minima takes place. We introduce notation uj 2 = (1 + 7)2 = 1 + 5. 

In order to calculate the amplitude (0\e~ TH \a) with instanton calculus, we start with formula analogical 



to (A.22): 



(a\e- TH \0) =e- SE ^Afdet- 



dr 2 



+ V"(x(r)/a) 



(A.92) 



where x(r) is a classical trajectory which starts at and ends at a. The potential V(x) satisfies conditions 
V"(0) = uj 2 and V"(l) = 1 which gives us information on asymptotic behavior of x(t) for infinite T. In 
analogy to formula (A. 36 1: 



x(t) = a — Ce *, t -> 
x{T) = C'e u \ t->-o 



(A.93) 
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We say that position of an instanton x{r) is t = tq if x(tq) = a/2. Let So = limr-yoo Se[x(t)] be 
approximate action of a single instanton. Due to the fact that action is symmetric under time rever- 
sal and parity transformation x — > —x, it is the same for instantons going between neighboring minima: 
(a,0), (0,a) (0,-a), (-a,0). It yields 



S = a z I dx^2V(x). 
la 



(A.94) 



The complete amplitude (a|e Tff |0) is a sum over many instanton solutions each of which has a form of 
RHS of (A.92I. A trajectory which begins at x = and ends at x = a jumps between various minima an 



odd number of times, so it contains and odd number of instantons. 



a e 



-TH 



|0)=5>: 



2n+l t 



-S B [2 2n+1 (T)]jy det - 



dr 2 



+ V"(x 2n+1 ( T )/a) 



(A.95) 



where x 2n+ i(r) is a classical trajectory which consists of 2n + l glued instantons. Positions of the instantons 
are t± < t 2 < . . . < T2 n +1- For r € (— CO, Ti) the trajectory is exponentially close to the middle minimum 
x = 0. Then it jumps to i = ±a and for r G (721+1, ^21+2) it is close to one of side minima. From the middle 
minimum it goes to the left or right minimum, so for r G {tu, T21+1) we have X2 n +i(T) ~ ±0. At last it goes 
to the right minimum so that X2u+i{t) ~ a for r G (t2„+i, 00). Because of the freedom of side minimum 
to which the instanton goes at r 2 i+i there is an additional coefficient A 2n +i = 2™ which is the number of 
possible paths. Positions of instantons correspond to zero modes of Euclidean action and they have to be 



integrated out. Normalization of zero modes gives an additional coefficient (\/5o/27r) 



2n+l 



He- Tff |0}=£2 



n e -(2n+l)S 



n 




2n+l 



< / dn ... dr 2n+ iAf [ det ' 



dr 



2 +V"(x 2n+1 (r)/a) 



-1/2 



(A.96) 



In analogy to formula (A. 50 1 we can write 

d 2 



N det 



dr 



2 +V"(x 2n+1 (T)/a) 



-1/2 



= (x = 0\U 



T 



\x = 0), 



(A.97) 



where U is the evolution operator 



I 'l t>. r, ! =- T «>xp ( jf 2 dr^P 2 + il/"(5 2 n+i(r)/a)X 2 



(A.98) 



Remember that the zero modes were integrated out, so we will be interested in calculating the expression 



N det ' 



dr 2 



+ V"(x 2n+1 (T)/a) 



-1/2 



= > = 0!r(|,-| ) k = 0}.V A,; (A.!)-)] 



where Ao is the zero mode corresponding to the instanton connecting the central minimum of V with left or 
right minimum and Ao to the instanton ending at the middle minimum. 

For r G (t2i+i + l,T2i+2 — 1) the trajectory x 2ti+ i(t) is in one of the side minima of the potential and 
V" ' (x 2n +i{ T ) 1 0) ~ 1 an d for t G (T 2 i+1, T2i+i— 1) it is in the central minimum at which V" {x 2u+ i{t) / a) ~ cu 2 . 
Thus, we introduce approximate evolution operators 



",.-,): J. xp ( (r / -r I )^P 2 + ^-X- 

/■ M-(^/-^)V-2* 2 



(A. 100) 



A.2. INEQUIVALENT MINIMA 



59 



and observe that Z7(r 2 i+i+l,r 2 i+2-l) ~ C/i(r 2l+ i+l, rai+a— 1) and U(T 2i + l, r 2 j+i-l) t/ w (r 2 i+l, r 2i+ i-l). 
Let us denote by \E^} and \E%) the eigenbases of harmonic oscillators with frequency 1 and uj respectively. 
We do now analogical transformation to ( |A.54 1 . In this case it is more complicated due to different harmonic 
frequencies in minima of the potential. 



T 



T 



= U r 2 „ +1 + lj . . . U (r 2 + 1, t 2 - 1) U (r 2 - l lTl + 1) U (n + l,7i - 1) £/ \n - 1, - ^ 
« ^> ra„ +1 + l) ■ • • f/ (r 2 + 1, 7 2 - 1) U x (r 2 - 1, r x + 1) U (n + 1, n - 1) [/„ (n - 1, -| 
= Ui (|,ra„ +1 + l) ■ • • E TO W ( T 2 + 1^2-1) E l^m) «l Ui (75 - l,n + 1) X 

x E TO TO ^ (n + i,n - i) E TO TOI ^ ( r i - !. ~f 

« Ui ^,r 2n+1 + l) TO . . . |£ W ) C/ (7 2 + 1,7 2 - 1) TO TO Ui(r 2 - l,n + l)x 
x l^ 1 ) [/(n + 1,7!- 1) \E%) (E%\ U u (t! - 1, -|) 

= Ul (|,7 2 „ +1 + l) TO • • • TO (E%\ ^ (7 2 + 1, 7 2 ) U X (7 2 , 7 2 - 1) |4) X 

x (^| u x (7 2 - i,7i + 1) TO {El\u x (t, + i,n) c/ w (n.n - 1) x 



x (sg'iCL n-i 



«(T2n+l) K/ ( r 2n) • • ■ k'(7 2 )k(7i) 



= e - ( ?-^ +l)£ ° |^> . . . \Eg) {E-\El) e-^-^ E o (Eq\Eq) (E*\E%) e-^ + ^ E o (E£\ x 

k(7 2 „ + i)k(7 2 „) . . . k(t 2 )k(ti), 



where 

In the last step we used relations 



(£ 1 |E/(7 + 1,7 -l)|£ff) 
(S 1 |C/ 1 (7 + 1,t )U u (t ,t - l)\E%) ' 
(grpl 1,7 -1)|^) 

(E$\U U (T + 1,70)^(70,70 - 1)1^) ' 



(A.101) 



(A. 102) 



(A. 103) 



U 1 {T l+1 ,r l )\El)=e-^-^ E °\El), 

The eigenbasis for the harmonic oscillator is known and one can easily calculate scalar products standing in 
the equation (A. 101 1 



(x = Ol^o 1 ) 



-1/4 



/UJ\ 1 / 4 

(rr = 0|K) = (-) , 
( J B 1 |K)= W 1 / 4 . 



(A. 104) 



1 + w 

Let us denote by T w the time which the trajectory spends in minimum with V"(x) = co and Ti the time 
which the trajectory spends in minimum with V"(x) — 1, i.e. left or right minimum, neglecting sizes of 
instantons: 



T u = T 2n+ i - r 2n + T 2n -i - r 2n _ 2 + . . . + ti - (-T/2), 

Ti = T/2 — T 2n+ l + T 2n — T 2n -1 + T 2n -2 + ... — 71. 



(A. 105) 
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Then there is T 1 +T LJ = T. Then, 



(x = 0\U 



T T 
2~'~~2 



\x = 0) 



s/tt 
<Jtt 



e -T„K e -TiEi 



e -T/2 e -lT„ 



,1/4 



1 + UJ 



2n+l 



t=0 



,1/4 



1 +0J 



2n+l 
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Recall that u? = (1 + 7) 2 and E% = uj/2. Terms k are corrections coming from the fact that a potential 
with varying frequency was approximated by a potential which is constant on intervals (t^, Ti + i). They can 
be written as 



k(to) 



det 



dr 2 



W(t) 



1/2 



det [-£ J + V"(x 1 (r)/a)} 
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where x\{t) is a classical trajectory consisting of one instanton located at tq going from x — to x — a. The 
function W(t) is the frequency which was used for approximation. Coefficient k corresponds to an instanton 
starting at middle minimum and ending at side minimum while k corresponds to minimum going from side 
to central minimum. For an instanton starting at the middle minimum we have 



According to (A.71I 



W(t) = 



U> T < T , 
1 T > T . 



TO 

MT/2) 



where 



+ V"(x 1 (T)/a))Mr) = 



dr 2 

M-T/2) = 0, 
M-T/2) = 1. 

with V" (xi(t) / a) replaced by W(t) for 4>o( T )- It is easy to show that 



isinh( W (T + T/2)) 
J sinh(w(T + T/2)) cosh(T - t ) + cosIi(w(t + T/2)) sinh(T - t ) t > t 



T < T 



(A.108) 

(A. 109) 

(A.110) 
(Alll) 

(A. 112) 



and therefore $}(T/2) « |(1 + k) e T+y(r +T/2) _ 

We will now find value of the function ipoi T ) a t point r = T/2. The one instanton trajectory has 
asymptotic behavior 



x(t) « aA + e 1 



i — > 00 
t — > —00 



(A.113) 
(A. 114) 



By translating argument of x(r) one can make the coefficients A + and A- equal: A + —> A' + = A and 

jo;/ 



A_ A'_ = A where A = ^/(i+<-) 41/(1+^ ne of solutions of dA.110) is j/ x (t) = -A-£(t). Its 



asymptotic behavior is 



t — > 00, 
i — > —00 



(A.115) 
(A.116) 
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where B = -7=-. We can normalize the other solution t/ 2 (t) of (A. 110 1 so that the Wronskian satisfies 

v So 



relation 

W = yizj 2 - yiJ/2 = 2B 2 . 
One can read off from Wronskian the asymptotic behavior of function 2/2(7"): 

2/2(7) -.Be*, t-toc, 



2/2 (t) « --Be 
10 



t — > —00. 



(A.117) 



(A.118) 



A function i^o(t) satisfying both, equation (A. 110 1 and boundary coditions (A. Ill I can be approximated by 
a specific linear combination of two solutions y\{r) and j/3(t): 



Then, 



^o(r) = ^(^ T/ V(t) + We - T / 2 2/ 2 (r)). 



v ' o(r/2) = i (e7T/2+we " 7T/2) 



(A.119) 
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From what was shown in Section [A. 1| we know that tp\ (T/2) = 0. Similarly as before, we have 
0«Vo(r/2)+ / 1 dr'G(T/2,T')X Mr') 

J T /2 

= MT/2) + / T/2 dr'^ 2 / 1 (r') 2 / 2 (T/2)A -l-(e" T / 2 y 1 (r') + We - T / 2 y 2 (r')) 



t/2 
T/2 



2cjB 



(A.121) 



The second approximation holds since e uT / 2 J rfr'?/i(T') 2 > e-" T / 2 J dT , y 1 (T')y 2 {T'). The last equality holds 
because of normalization of y\ (r) . We conclude that 

A = 4ujB 2 e- T -^ T il; (T/2) (A. 122) 

We are interested in calculating the quantity 



k(t )V^o 



/Vo (T/2)A 



Be 7To/2 vTTw 



/ ^o(T/2) 

One can calculate k(tq)\/Xo precisely in the same manner. The difference is that now 

W(t) = 

so the functions 2/1 (r) and y 2 are different than previously. The result yields 

«(to)\/^o = Be~ 1To/2 V%j 
Finally, we can calculate right side of ( |A.99 1 



1 t < T 

LU T > Tq 



(A. 123) 



(A.124) 



(A. 125) 



Af (det '[...])" 



1/2 = ^ e -T/2 e - } T„ fl/4 



2n+l 



1 + U) J 

n n 

x]jBe^ /2 VTT^]jBe-^+ l/2 V2^ 

i=0 i=l 

^(T^f V * <,+i, ( 23/4B ^< 1 + ^ 



(A. 126) 
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Applying (A. 96) we get the final formula 



a e 



-TH 



|0) ~S (2n + l)!V27r (l + w 



1/4 



e -|(i+i) 2 3/4 



w5 



2n+l 



-So 



(l+w)- 1 / 4 BT 



27T V 1 + CJ 



Recall that 5 = where A 

vS 



1/4 



: ( 1+ *)sinh ^y^e-^Cl+wJ-^flrj 



(A. 127) 
(A. 128) 



= ^7(1+^)^1/(1+^) 



two energies: 



and A± are defined in (A. 113 1. This amplitude reveals 



£-^(l + ^2 3/4 (l + .r 1/4 ^eAA 
£ + = \{l + \) + 2 3/4 (l + u)-^^e- s °aA. 



(A. 129) 



-1/2 



Let us now discuss general expression for the determinant 7V( det ' [ — + V"(x„(r) /a) J ) " when the 
classical trajectory x„(t) connects two minima with V"(x n (— oo)) = Wi and V"(xn(+oo)) = W2- For wi = W2 
the number of instantons which build the trajectory x n (r) must be even, i.e. n is even. For wi ^ W2 number 
of instantons n is odd. For each instanton going from central minimum to one of side minima there is the 
lowest eigenvalue corresponding to zero mode %/Ao an d the term k(tj). For each instanton going the other 
way there are similar terms \/Ao and k(tj). For any transition there is additionally (Eq\Eq). The Euclidean 
evolution of a state in the middle minimum gives e~ TujE ° where T u is the total time which the trajectory 
spends in the central minimum. Evolution in side minima gives e~ TlE ° where T\ = T — T w . Similarly as 
before there is also term (x = 0\Eq 2 ) (x = O^ 1 ) = (u^) 1 / 4 /^. If uji = UJ2 then there are n/2 instantons 
going from cental to side minimum and n/2 instantons going from side to central minimum. If lo\ = 1 and 
uj2 — w there are (n+l)/2 instantons starting at side minimum and ending at central minimum and (n— 1)/2 
the other instantons. Inserting expressions for «;(rj) and Aq we get 



AM det ' 



'dr 2 



+ V"(x n (r)/a) 



(0J1UJ2) 



1/4 



,1/4 



aA 



+ UJy/So I 

(2w(l+a;))"/ 4 

(2w(l+a;))(™- 1 )/ 4 yi _ KJ 

(2w(l+w))("- 1 )/ 4 y2^ 



w a = a; 2 

Ul\ = U), Ll>2 = 1 
Wi = 1,^2 = W 
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It is now simple to calculate other amplitudes, e.g. (a\e~ \a). The difference is that now there are even 
number of instantons. An 2n instanton solution has 2"~ 1 topologically different trajectories so A^n = 2 ,l_1 . 
The only exception is N — 1 because there exists a trivial constant trajectory. We may write N 2n — 



+ id n Q. According to (A. 130 1, with u)\ = = 1) 



A/^det' 



dr 2 



+ V"(x 2n (T)/a) 



-1/2 



2 3/4 ^v^(i+-r 1/4 ) 



2)i 



In analogy to (A.96) there is 



^ (2n)! I V 2tt 




dr\ . . . dT2 n N ( det ' 



+ V"(x 2n {T)/a) 



1 



2^r 



e 2 cosh 
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There are three energies: 



1 2 

E 2 = 1 -+2^(l+oj)- 1 /\l^e- s °aA 



(A. 133) 



We can see that this result is different than (A. 129 1. Still, this can be explained. The first thing is that 



when we calculated (0|e TH \a) there were two energies. It is because (Ei\0) = and the expansion of the 



amplitude 



-TH 



a) is 



\e- TH \a) 



J2^ TE "(0\E n )(E n \0) 



(A. 134) 



Secondly, energies calculated using different amplitudes differ by a perturbative quantity J! We know 
however that the instanton calculus does not take into account perturbative corrections of energies and thus 
any perturbative terms appearing in semiclassical approximation in Euclidean space can shall be ignored. 
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